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ABSTRACT 


A  wireless  sensor  network  (WSN)  is  eomprised  of  sensor  nodes  designed  to  collect  and 
transmit  data  efficiently.  For  this  reason,  WSNs  are  relied  upon  by  the  Department  of 
Defense  for  deployment  in  remote  and  hostile  areas.  The  performance  of  a  WSN  is 
degraded  by  the  amount  of  interference  experienced  by  nodes  during  simultaneous 
transmissions.  Transmitting  in  the  presence  of  interference  can  affect  the  lifetime  of 
sensor  nodes  by  requiring  multiple  re-transmissions  of  data.  In  this  thesis,  we  propose  a 
routing  algorithm  that  uses  spatial  time-division  multiple  access  (STDMA)  to  schedule 
simultaneous  transmissions  such  that  interference  is  mitigated  and  transmission  time  slots 
are  reused  appropriately.  We  integrate  STDMA  with  a  physical  interference  model  that 
facilitates  the  calculation  of  interference  metrics  based  on  signal-to-interference  ratio. 
Using  the  interference  metrics  as  link  costs,  we  implement  Dijkstra’s  algorithm  to 
determine  the  least  interference  path  from  a  sensor  node  to  the  gateway.  Via  simulations 
using  MATLAB  and  QualNet,  we  show  that  this  approach  to  interference  mitigation 
helps  network  performance  by  decreasing  end-to-end  delay.  We  develop  this  algorithm 
as  a  proof-of-concept  to  show  that,  despite  the  computational  complexity  associated  with 
interference  based  scheduling,  STDMA  can  have  a  real  impact  on  network  design  and 
performance. 
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EXECUTIVE  SUMMARY 


Efficient  design  and  implementation  of  wireless  sensor  networks  (WSNs)  has 
become  an  increasingly  active  area  of  research  in  recent  years  due  to  the  vast  potential  of 
sensor  networks  to  enable  applications  that  connect  the  physical  world  to  the  virtual 
world  [1].  There  is  signifieant  interest  in  the  development  and  use  of  WSNs  for  military 
applications  to  increase  battlefield  information  dominance.  WSNs  provide  the  ability  to 
maintain  an  information  gathering  eapability  in  remote  and  hostile  environments. 

Like  all  wireless  networks,  the  performance  of  a  WSN  is  degraded  by  the  wireless 
interference  experieneed  at  individual  nodes  during  simultaneous  data  transmissions. 
Additionally,  WSNs  have  limitations  due  to  the  limited  battery  capacity  of  the  nodes  [1]. 
The  lifetime  of  a  sensor  node  can  be  signifieantly  reduced  if  exeessive  wireless 
interference  neeessitates  multiple  re-transmissions  of  data.  Thus,  interferenee  and  sensor 
node  lifetime  go  hand  in  hand.  When  designing  a  network  protocol  stack  suitable  for  use 
in  WSNs,  interference  at  the  physical  layer  must  be  taken  into  aceount  to  provide  for 
reliable  and  effieient  WSN  operation. 

In  this  thesis  we  eoncentrate  on  modeling  interferenee  by  exploiting  a  eoncept 
known  as  spatial  reuse.  Essentially,  spatial  reuse  allows  nodes  to  simultaneously  transmit 
during  the  same  time  slot  as  long  as  they  have  sufficient  spatial  separation.  Spatial  reuse 
is  modeled  at  the  data  link  layer  (DLL)  using  spatial  time-division  multiple  aeeess 
(STDMA)  [2].  STDMA  governs  the  sehedule  of  when  nodes  transmit  (i.e.,  which  time 
slot  is  assigned  to  a  node  for  transmission).  When  using  STDMA,  multiple  nodes  are 
assigned  to  transmit  during  the  same  time  slot,  thus  allowing  appropriate  time  slot  reuse 
to  increase  network  performance.  We  integrate  the  STDMA  protoeol  with  a  physical 
interference  model  in  order  to  capture  interferenee  metrics  using  signal-to-interference 
ratio  (SIR).  Our  focus  is  to  use  these  captured  interference  metrics  to  determine  the  least 
interfering  paths  between  sensor  nodes  and  the  gateway  to  allow  for  mitigated 
interference  routing. 
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In  order  to  solve  the  problem  of  routing  in  a  WSN  with  interference  generated 
from  simultaneous  transmissions  (inter-node  interference),  it  is  important  to  first 
understand  the  models  that  exist  to  capture  interference  in  wireless  networks.  The  nature 
of  the  wireless  communications  medium  used  in  WSNs  results  in  the  inclusion  of  inter¬ 
node  interference  in  addition  to  the  desired  transmitted  signal  at  all  receivers  in  the 
network.  Typically  SIR  is  the  key  parameter  in  evaluating  signal  reception  in  wireless 
networks.  There  are  two  primary  interference  models  used  in  networking  literature  to 
capture  the  effects  of  interference  on  wireless  transmissions.  These  are  the  Protocol 
Interference  Model  (PrIM)  and  the  Physical  Interference  Model  (PhIM).  Under  the  PrIM, 
two  nodes  can  successfully  communicate  if  and  only  if  the  receiver  is  within  the 
transmission  range  of  the  sender  and  the  receiver  is  out  of  the  interference  range  of  other 
simultaneous  senders.  Under  the  PhIM,  a  receiver  can  successfully  receive  data  from  the 
sender  if  and  only  if  the  SIR  of  the  sender  at  the  receiver  is  no  less  than  a  predefined 
threshold.  The  PrIM  simplifies  the  communication  model  of  wireless  networks  and  is 
therefore  more  convenient  for  analysis.  By  contrast,  the  PhIM  takes  the  received  physical 
signal  strength  as  a  criterion  to  measure  interference,  which  has  been  proven  to  be  more 
accurate  and  reliable  [2]. 

In  this  thesis,  we  use  the  PhIM  to  capture  interference  metrics  in  the  WSN. 
Specifically,  we  use  STDMA  to  schedule  transmissions  and  then  use  PhIM  to  calculate 
the  SIR  witnessed  at  each  receiving  node. 

We  implement  STDMA  with  the  following  in  mind:  Given  our  goal  of 
maximizing  spatial  reuse  and  controlling  interference,  we  do  not  allow  one -hop  neighbors 
to  simultaneously  transmit  because  of  the  excess  interference  that  would  be  created  from 
these  transmissions.  We  allow  for  two-hop  neighbors  to  transmit.  These  neighbors  form 
a  transmission  group  (i.e.,  all  nodes  within  the  same  transmission  group  can  transmit 
simultaneously).  An  example  network  with  color-coded  transmission  groups  is  shown  in 
Figure  1.  For  example,  in  Figure  1,  all  blue  nodes  may  transmit  simultaneously  in  the 
same  time  slot.  This  holds  for  all  the  colored  node  groups. 

The  next  step  is  to  compute  the  interference  experienced  during  each  reception 
due  to  simultaneous  transmissions  throughout  the  transmission  group.  In  other  words. 


what  interference  does  node  2  experience  when  listening  to  node  1  due  to  transmissions 
from  nodes  3,  5,  11,  13,  and  15?  For  this  calculation  we  assume  that  all  nodes  in  the 
transmission  group  have  data  to  send,  which  provides  the  maximum  possible 
interference.  The  calculations  for  the  SIR  at  each  node  can  be  calculated  as 


SIR,U)  = 


P.iJ) 

S  PM ' 

k^TGj 


(1) 


Network  topology 


Figure  1 .  Example  of  a  network  that  depicts  color-coded  transmission  groups.  All 
nodes  within  the  same  color-coded  transmission  group  may  transmit  simultaneously 

within  the  STDMA  schedule. 


The  SIR  calculation  given  by  Equation  (1)  states  that  the  interference  received  at  node  i 
during  reception  from  transmitting  node  j  (i.e.,  node  2  from  node  1)  is  the  aggregate  of 
the  transmitted  power  from  all  other  nodes  in  the  transmission  group  of  node  j  (TGj), 
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which  is  received  at  node  i .  Equation  (1)  maintains  the  assumption  that  all  nodes  in  the 
transmission  group  send  data  simultaneously  and  that  the  interference  is  at  its  maximum 
possible  value.  This  caleulation  is  repeated  for  every  eombination  of  reeeiving  nodes  and 
applieable  transmission  groups,  which  enables  us  to  assign  an  interference  metric  (i.e., 
the  caleulated  SIR  value)  to  each  edge  in  our  interference  graph. 

To  faeilitate  the  SIR  caleulation,  we  take  advantage  of  graph  theoretie  models. 
Any  network  ean  be  represented  as  a  graph  G(V,  E)  where  V  is  the  set  of  nodes/vertiees 
in  the  network  and  E  is  the  set  of  edges/links  in  the  network.  We  transpose  the  original 
network  graph  G  into  what  is  known  as  a  conflict  graph  Gi.  In  a  conflict  graph,  denoted 
asG^(F/,E'/) ,  a  vertex  is  introduced  for  eaeh  link  in  the  original  network G(E, £')  .  An 

edge  in  Gj{V^,Ej)  eonnects  two  vertiees  in  the  conflict  graph  if  these  two  links  interfere 
in  the  original  graph  (as  in  the  ease  of  simultaneous  transmissions).  To  characterize  the 
interference  relationship  between  two  links  in  the  eonflict  graph,  we  assign  a  SIR  to  the 
links  that  refleet  the  aggregate  SIR  for  all  simultaneously  transmitting  nodes.  We  treat 
the  SIR  as  link  weights  such  that  a  higher  SIR  equates  to  a  lower  cost  for  transmission. 
Specifically,  we  assign  to  each  link  in  the  conflict  graph  the  value  of  the  inverse  SIR  for 
the  receiving  node  of  the  link.  The  graphical  results  of  the  link  weight  calculations  and 
assignments  for  the  network  in  Figure  1  are  shown  in  Figure  2. 

We  use  the  link  weights  and  connectivity  graphs  developed  above  as  inputs  to 
Dijkstra’s  least  cost  algorithm  to  calculate  the  shortest  path  from  a  sensor  to  the  gateway 
with  least  interferenee. 

Using  MATFAB  and  QualNet,  we  found  that  the  algorithm  we  developed 
improves  network  delay  but  hinders  network  throughput  for  various  fading  and  traffie 
models  when  eompared  to  the  traditional  time-division  multiple  access  (TDMA)  scheme. 
The  results  shown  in  Table  1  compare  the  network  performance  for  several  scenarios  that 
use  a  TDMA  sehedule  and  several  scenarios  that  use  interference-aware  STDMA 
scheduling.  In  addition  to  different  transmission  sehedules,  the  scenarios  include  variable 
packet  traffic  types  and  fading  ehannels  for  both  TDMA  and  STDMA  seheduling.  The 
results  of  Table  1  specifieally  relate  the  performance  of  TDMA  versus  STDMA  in  terms 
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of  throughput  and  delay.  From  the  results,  it  ean  be  seen  that  there  is  a  signiHeant 
improvement  to  network  delay  when  using  STDMA  rather  than  TDMA.  However,  the 
throughput  performance  does  not  increase  notably  and  in  some  configurations  decreases 
or  remains  the  same. 

Our  interpretation  of  the  results  is  that  time  slot  re-use  allows  slots  to  matriculate 
faster  across  the  network  but  at  the  expense  of  interference  causing  the  received  SIR  to 
drop  below  the  required  threshold.  This  loss  of  packets  results  in  a  decreased  throughput. 

Despite  the  tradeoff  between  delay  and  throughput,  we  see  the  results  obtained  as  a 
proof-of-concept  that  a  deterministic  PhIM  can  provide  value  in  network  design  and 
performance  despite  its  computational  complexity. 

An  interference  aware  routing  mechanism  that  considers  the  scheduling  of 
simultaneous  transmissions  using  STDMA  can  offer  quantifiable  gains  to  the  overall 
network  performance  of  WSNs.  The  technique  presented  in  this  thesis  as  a  proof-of- 
concept  shows  that  the  PhIM  integrated  with  a  STDMA  approach  to  obtain  least 
interfering  paths  can  provide  for  efficient  utilization  of  network  resources  for  a  given 
network  implementation. 


Table  1 .  Results  comparing  delay  and  throughput  when  using  STDMA 
versus  TDMA  under  different  fading  and  traffic  models  for  a  10-node  network. 


MAC 

Protocol 

Fading 

Model 

Traffic 

Model 

Throughput 

(kbps) 

Delay 

(sec) 

Throughput 
Improvement 
using  STDMA 
(%) 

Delay 

Improvement 
using  STDMA 
(%) 

TDMA 

None 

CBR 

4283 

0.1237 

-0.4 

24.3 

STDMA 

None 

CBR 

4266 

0.0937 

TDMA 

Rayleigh 

CBR 

4105 

0.1194 

3.9 

21.5 
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Figure  2.  Conflict  graph  with  associated  SIR  values  for  the  network  graph 
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XXIV 


I.  INTRODUCTION  AND  MOTIVATION 


In  recent  years  there  has  been  a  rising  interest  in  data-centric  warfare  in  all 
branches  of  the  military.  Wireless  sensor  networks  (WSNs)  provide  a  flexible 
technology  option  to  achieve  the  objective  of  gathering  information  about  the  physical 
environment  and  leveraging  this  information  to  achieve  battlespace  dominance  in  the  air, 
on  land  and  at  sea.  The  miniaturization  of  computer  processors  and  advancement  of 
communications  technologies  allows  for  a  widespread  distribution  of  sensors,  in  an  ad- 
hoc  manner,  for  data  extraction  in  real  time. 

WSNs  are  wireless  ad-hoc  networks  that  are  characterized  by  low  power,  low 
processing  capable  measurement  nodes  that  gather  data  from  the  environment.  The 
sensor  nodes  also  act  as  routing  elements  for  the  movement  of  this  data  from  a  source 
node  to  a  destination.  The  destination  node  in  a  WSN  is  characterized  as  a  gateway  or 
sink  node.  The  gateway  is  a  point  of  aggregation  where  information  from  all  the  sensors 
is  received  and  processed.  The  gateway  eventually  sends  this  gathered  information  to  a 
command  and  control  unit  for  processing  and  dissemination. 

The  small  physical  size  and  inherently  limited  capabilities  of  the  WSN  nodes 
combined  with  the  desire  to  use  the  nodes  without  replacement  serves  as  the  motivation 
to  develop  simple  and  efficient  communications  mechanisms  for  a  WSN.  The  most 
lifetime-limiting  component  for  a  WSN  node  is  the  power  supply,  and  the  most 
demanding  power  consuming  component  is  the  transceiver  antenna  in  the  transmitting 
mode.  It  is  with  these  considerations  in  mind  that  it  becomes  readily  apparent  that  an 
effort  must  be  made  to  minimize  transmissions  from  the  sensor  node  while  still 
forwarding  the  required  collected  information. 

In  a  WSN,  which  can  consist  of  hundreds  or  thousands  of  nodes,  sequencing  the 
transmission  and  movement  of  data  is  required  to  minimize  delay  of  data  delivery  from 
the  network  and  to  prevent  node-to-node  interference  as  data  flows.  This  unique  routing 
demand,  coupled  with  limited  node  battery  life,  requires  WSN  communications  protocols 
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to  be  tailored  to  provide  low  power  eonsumption  and  eongestion  free  data  transfer.  With 
these  protoeols,  WSNs  aeross  a  potentially  large  physieal  area  in  an  effort  to  provide  high 
end-to-end  data  throughput. 

This  thesis  foeuses  on  modeling  inter-node  interferenee,  and  using  the 
interferenee  relationships  as  a  basis  for  a  least-eost  routing  algorithm  to  inerease  the  data 
throughput  of  a  WSN.  For  the  remainder  of  this  thesis,  the  term  interferenee  refers  to 
inter-node  interferenee  that  is  generated  by  simultaneous  transmissions. 

A,  WIRELESS  SENSOR  NETWORK  NODE  ARCHITECTURE 

The  desire  for  aeeurate  and  timely  remote  gathering  of  information  has  served  as 
the  motivation  for  the  development  and  rapid  maturation  of  WSN  teehnologies.  WSNs 
serve  an  important  purpose  in  military  speeifie  applieations  sueh  as  foree  proteetion 
systems  designed  to  deteet  and  report  deteetion  of  enemy  movements. 

The  WSN  nodes  eonsist  of  four  key  eomponents  in  addition  to  any  number  of 
optional  applieation  dependent  eomponents.  The  key  eomponents  are  a  sensing  unit,  on¬ 
board  proeessor,  radio  transeeiver,  and  power  supply.  These  WSN  nodes  are  used  for 
environmental  sensing,  data  proeessing,  communieations,  and  stand-alone  battery  power 
[1].  The  key  eomponents  of  a  WSN  node  are  shown  in  Figure  1.  The  sensor  on  a  WSN 
node  is  designed  to  sense  a  speeifie  physieal  quantity,  transfer  this  physieal  quantity  into 
an  eleetrieal  signal,  and  then  digitize  the  information  for  subsequent  proeessing  and 
forwarding.  Onee  this  digital  information  has  been  eolleeted,  the  node’s  proeessing  unit 
performs  a  variety  of  tasks  ineluding  data  averaging  and  paeket  eonstruetion  in 
preparation  for  transmission.  The  radio  transeeiver  transmits  paekets  that  originate  from 
the  node  while  also  forwarding  paekets  reeeived  from  surrounding  nodes.  The  power 
supply,  a  battery,  is  used  to  power  all  the  aforementioned  funetions  and  is,  therefore,  a 
eritieal  design  feature  of  any  WSN  node. 
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Figure  1.  WSN  node  arehiteeture  from  [1], 


B,  LAYERING  ARCHITECTURE  OF  A  NETWORK 

When  designing  a  eommunieations  network  for  widespread  interoperability,  the 
most  eommon  model  used  is  the  Open  Systems  Interconnection  (OSI)  reference  model 
shown  in  Figure  2.  The  OSI  model  is  used  a  reference  because  it  reduces  the  complexity 
of  network  design  and  analysis  and  separates  the  functions  logically  into  modular 
sections  to  create  an  understood  framework  for  network  implementation.  The  OSI  model 
consists  of  seven  distinct  protocol  layers  that  are  often  referred  to  as  a  protocol  stack. 
Each  layer  in  the  stack  communicates  with  the  layers  above  and  below  it.  The  layers  of  a 
network  operate  to  provide  reliable  and  effective  data  communication  over  a  variety  of 
physical  media  channels  based  on  the  needs  of  the  network  designer  and  the  network 
users.  The  lower  layers  of  the  model — physical  layer  (PHY),  data  link  layer  (DLL), 
network  layer,  and  transport  layer — are  primarily  concerned  with  the  formatting, 
encoding,  and  transmission  of  data  over  the  network  [2].  The  functions  of  these  layers 
are  implemented  in  both  hardware  and  software  in  WSNs.  The  higher  layers  of  the  OSI 
model — session,  presentation,  and  application — are  primarily  concerned  with  interacting 
with  the  user  and  are  usually  implemented  in  software  [2]. 
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Figure  2.  OSI  network  layering  model  from  [3]. 


The  typical  WSN  networking  protocol  stack  consists  of  a  number  of  layered 
protocols  similar  to  the  OSI  network  stack.  Specifically,  there  is  a  physical  layer  (PHY 
Layer)  for  radio  communication,  data  link  layer  (DLL)  for  packet  transmission 
coordination,  and  a  networking  layer,  which  provides  an  efficient  network-wide  packet 
routing  solution. 

The  PHY  layer  specifies  all  network  communications  parameters  such  as 
frequency  bands  and  modulation  schemes.  In  order  for  the  PHY  layer  to  operate 
correctly,  each  node  must  use  the  same  pulse  shaping  and  timing  mechanisms.  Choices 
for  PHY  layer  parameters  are  based  on  balancing  the  needs  of  providing  adequate 
performance  while  limiting  the  complexity  and  power  required  for  operation. 

The  DLL  is  typically  separated  into  the  medium  access  control  (MAC)  and  logical 
link  control  (LLC)  sub-layers.  The  MAC  sub-layer  serves  to  allow  the  users  coordinated 
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access  to  the  transmission  medium  by  the  use  of  either  distributed  or  eentralized 
eoordination  funetions.  Distributed  eoordination  involves  the  use  of  information  loeal  to 
eaeh  node.  Nodes  then  make  logical  decisions  in  order  to  aehieve  eollision  free 
transmissions.  Centralized  eoordination  deseribes  the  use  of  a  eentralized  entity  to 
eoordinate  aceess  to  the  transmission  medium.  In  the  ease  of  a  WSN,  the  entity  is  the 
gateway  node  and  is  used  to  ereate  a  transmission  sehedule  that  is  sent  to  all  other  nodes. 
An  example  of  a  MAC  seheduling  funetion  is  the  time-division  multiple  aeeess  (TDMA) 
seheme  wherein  eaeh  node  is  assigned  a  time  slot  on  a  sehedule  during  whieh  it  is 
allowed  to  transmit  information  to  neighboring  nodes.  TDMA  allows  eollision  free 
transmission  with  little  eomplexity  but  does  not  maximize  network  performanee  beeause 
only  one  node  ean  transmit  at  a  time,  and  eaeh  node  remains  idle  for  long  periods  of  time. 
In  order  to  improve  the  performanee  of  TDMA  seheduling,  spatial  time-division  multiple 
aceess  (STDMA)  allows  multiple  nodes  to  utilize  the  same  time  slot  if  they  are  loeated  at 
a  suffieient  distanee  away  from  eaeh  other  in  order  to  limit  the  interference  eaused  by 
simultaneous  transmissions. 

The  network  layer  is  used  to  establish  the  optimal  end-to-end  routing  path  aeross 
the  network.  Additionally,  the  networking  layer  provides  eomponent  (i.e.,  node  and 
gateway)  addressing,  packet  fragmentation,  and  packet  reassembly. 

There  exists  a  wide  variety  of  Institute  of  Eleetrieal  and  Eleetronies  Engineers 
(IEEE)  standards  that  ean  be  used  as  the  design  basis  for  a  reliable  wireless 
eommunieations  system.  In  Chapter  II,  we  deseribe  the  IEEE  802.15.4  standard,  whieh 
speeifies  the  PHY  and  DEE  layers  for  low-rate  personal  area  networks  (ER-PAN).  This 
standard  provides  the  lower  layer  implementation  definitions  for  many  WSN  protoeol 
designs  and  is,  therefore,  an  appropriate  example  to  provide  a  framework  for  the 
funetionality  of  the  PHY  and  DEE  layers  in  a  WSN.  There  are  a  number  of  network 
layer  protoeols  that  are  appropriate  for  WSN  operation.  We  ehoose  the  Dijkstra  least- 
eost  routing  algorithm,  whieh  is  discussed  in  Chapter  V. 

In  addition  to  the  three  lower  layers,  there  are  protoeols  layered  above  these  that 
provide  applieation  speeifie  funetionality.  In  this  thesis  we  are  not  eoneemed  with  the 
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functionality  of  these  top  layers.  Our  focus  is  on  the  design  and  implementation  of 
routing  protoeols  considering  the  PHY,  DLL,  and  network  layers. 


C.  INTERFERENCE  MODELING 

As  wireless  communication  systems  have  evolved,  with  a  more  aggressive 
utilization  of  spectral  resources  and  more  effieient  transmission  systems,  a  better 
understanding  of  the  effects  of  interference  on  the  performance  of  sueh  systems  is  of 
paramount  importance.  This  has  motivated  the  development  of  more  accurate 
interference  models  [4].  The  resource  constraints  of  WSN  nodes  serve  as  a  motivating 
factor  to  design  the  most  efficient  communieations  protoeols  at  the  PHY  layer  in  an  effort 
to  increase  performance  at  higher  layers  in  the  presence  of  interference. 

The  purpose  of  an  interference  model  is  to  determine  the  effects  of  interference  at 
different  protocol  layers.  An  interference  model  can  be  viewed  as  a  combination  of  the 
following  sub-models: 

•  Propagation  Model:  Uses  eommunications  theory  to  deseribe  the  effects 
of  attenuation,  shadowing,  noise,  and  fading  on  the  received  signal. 

•  Transmitter  Distribution  Model:  Describes  how  the  transmitters  (i.e., 
nodes  in  a  WSN)  are  distributed  in  the  area  of  interest. 

•  Network  Operation  Model:  Deseribes  the  MAC  access  teehnique  in  the 
network,  whieh  specifies  the  network  eoordination  function  (i.e.,  the  use 
of  a  transmission  sehedule). 

•  Traffic  Model:  Describes  the  assumed  statistical  distribution  of  arrival  of 
packets  into  the  network. 

The  components  of  the  sub-models  appear  in  all  interference  models  and  can  be 
modeled  as  deterministic  or  random  processes  according  to  the  scenario  considered. 
Interference  models  are  largely  divided  into  (1)  statistical  interference  models,  which 
attempt  to  statistically  describe  the  wireless  interference,  and  (2)  interference  models  that 
seek  to  describe  the  effects  of  the  interference  on  the  operation  of  the  network  of  interest 
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[4],  The  two  most  prominent  subdivisions  of  the  latter  group  are  the  protoeol 
interferenee  model  (PrIM)  and  physieal  interferenee  model  (PhIM).  In  addition  to  the 
PrIM  and  PhIM  models,  there  are  graph-based  models  that  have  been  developed  to 
deseribe  the  effeets  of  interference  in  the  context  of  transmission  scheduling  and  node 
topology  control.  The  PrIM  states  that  a  transmission  between  two  nodes  is  successfully 
received  if  the  receiving  node  is  within  the  transmission  radius  of  the  transmitting  node 
and  outside  the  interference  radii  of  all  simultaneously  transmitting  nodes.  The  PhIM 
states  that  transmission  between  two  nodes  is  successfully  received  if  the  signal-to- 
interference  ratio  (SIR)  is  greater  than  a  minimum  value  [5].  In  general,  the  PrIM  allows 
for  faster  interference  computation  and  the  PhIM  provides  more  accurate  interference 
modeling.  All  classes  of  interference  models  relate  to  the  radio  capture  phenomenon, 
which  describes  the  ability  of  a  node  to  successfully  receive  a  radio  transmission  in  the 
presence  of  inter-node  interference  [4]. 

D,  MOTIVATIONS  AND  CONTRIBUTIONS  OF  THE  THESIS 

In  the  wireless  networking  literature,  there  are  many  studies  that  validate  the 
advantages  of  using  a  PhIM  approach  in  contrast  to  a  simpler  PrIM  for  interference 
modeling.  Despite  the  fact  that  the  PhIM  provides  a  more  accurate  representation  of 
interference  in  a  WSN,  this  accuracy  comes  at  the  cost  of  a  significant  increase  in 
complexity  for  the  MAC  and  PHY  layer  protocols.  In  this  thesis  we  propose  a  simplistic 
model  of  a  PhIM  in  order  to:  (I)  perform  a  weighted  conflict  graph  analysis  based  on  the 
SIR  on  all  links  in  a  WSN,  (2)  use  the  SIR  obtained  in  the  graph  analysis  as  the  cost 
metric  for  a  network  layer  routing  algorithm,  and  (3)  demonstrate  the  advantages  of  using 
a  least-cost  routing  algorithm  in  a  STDMA  MAC  protocol  in  comparison  to  a  TDMA 
MAC  protocol. 

The  contribution  of  the  research  proposed  in  this  thesis  is  a  proof-of-concept  for  a 
cross-layer  routing  and  medium  access  system  using  physical  interference  as  the  routing 
criteria.  The  WSN  models  and  assumptions  aid  in  the  demonstration  of  the  proposed 
concept.  Future  work  is  proposed  that  will  provide  a  more  robust  protocol  level 
implementation  for  a  more  general  WSN  setup. 
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E, 


ORGANIZATION  OF  THE  THESIS 


The  remainder  of  this  thesis  is  organized  as  follows.  An  in-depth  explanation  of 
the  OSI  reference  model  and  a  detailed  discussion  of  the  IEEE  802.15.4  Standard  for  ER- 
PAN  is  provided  in  Chapter  II.  The  concepts  of  interference  modeling  are  detailed  in 
Chapter  III  by  describing,  comparing,  and  contrasting  the  most  widely  used  interference 
models  in  the  literature.  The  MAC  protocol  implementation  and  the  advantages  of  using 
STDMA  in  lieu  of  TDMA  for  WSNs  are  described  in  Chapter  IV.  The  concept  of  least- 
cost  routing  and  the  operation  of  the  Dijkstra  algorithm  as  an  implementation  of  a  least- 
cost  routing  mechanism  is  explained  in  Chapter  V.  The  experimental  setup  and  an 
evaluation  of  the  simulation  results  is  discussed  in  Chapter  VI.  The  conclusion  and 
recommended  future  work  are  presented  in  Chapter  VII.  All  the  MATEAB  code  used  in 
the  simulations  are  contained  in  the  appendix. 
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II.  NETWORK  LAYER  PROTOCOLS 


In  this  chapter  we  present  an  explanation  of  both  the  general  eonstruetion  of 
layered  network  protoeols  and  the  speeifie  example  of  the  implementation  of  the  IEEE 
802.15.4  standard  for  WSN  implementation.  As  briefly  discussed  in  Chapter  I,  Section 
B,  modular  network  layering  protoeols  were  introdueed  by  the  OSI  referenee  model  in 
the  late  1970s.  This  approaeh  to  network  design  remains  the  standard  approaeh  to 
developing  a  network  protoeol  staek.  The  framework  provided  by  the  OSI  model  is  used 
in  nearly  all  WSN  designs  and  is  used  as  a  basis  for  discussion  of  the  network  layers  in 
this  seetion.  In  addition  to  the  general  information  provided  by  the  OSI  model  about  the 
fimetionality  provided  by  a  given  protoeol  layer,  the  IEEE  802.15.4  standard  is  used  as  an 
example  to  provide  a  WSN  specific  implementation  of  the  two  lowest  layers  of  the  staek. 
Beyond  the  PHY  and  DEE  layers,  the  discussion  of  higher  layer  functionality  in  this 
thesis  is  focused  primarily  on  the  routing  protoeol  of  the  network  layer. 

A,  PHYSICAL  LAYER 

The  lowest  layer  in  the  OSI  model,  the  PHY  layer,  is  responsible  for  the 
transmission  and  reeeption  of  data.  It  dietates  the  radio  frequency  bands  to  be  employed 
and  the  type  of  frequeney  spreading  and  modulation  techniques  to  be  used  for 
eommunication  [6].  The  original  IEEE  802.15.4  was  released  in  2003  and  has  received 
five  updates,  ineluding  the  most  reeent  update  in  201 1.  With  eaeh  of  these  releases,  the 
PHY  layer  has  been  updated  to  provide  more  flexibility  and  inereased  performanee  with 
an  aeeompanying  inerease  in  eomplexity. 

The  original  U.S.  IEEE  802.15.4  standard  provides  frequencies  in  the  2.4  GHz 
industrial  seientifie  medieal  (ISM)  radio  band,  for  unlieensed  deviee  operation,  using 
direet  sequenee  spread  speetrum  teehniques.  This  wide-band  implementation  supports  16 
ehannels  with  a  guard  band  of  5  MHz  between  them.  The  bands  range  from  2400-2483.5 
MHz  and  support  a  data  rate  of  250  kb/s.  The  seheme  uses  offset  quadrature  phase-shift 
keying  modulation  and  half-sine  pulse  shaping.  Updates  to  the  standard  have  inereased 
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the  bandwidth  from  wide-band  to  ultra  wide-band  (UWB),  which  increased  the  data  rate 
to  up  to  1  Mb/s.  Additionally,  multiple  modulation  schemes  were  added  to  allow  for 
variable  data  rates,  which  resulted  in  a  trade-off  between  data  transfer  rate  and  power 
consumption  [6]. 

B,  DATA  LINK  LAYER 

The  second  lowest  layer  in  the  OSI  model,  the  DLL,  is  primarily  responsible  for 
the  multiplexing  of  data  streams,  medium  access  control,  error  control,  and  data  frame 
detection.  As  was  stated  in  Chapter  I,  the  DLL  is  divided  into  two  sub-layers — a  higher 
LLC  sub-layer  and  lower  MAC  sub-layer.  The  subdivision  of  the  DLL  into  two  sub¬ 
layers  is  necessary  to  accommodate  the  logic  required  to  manage  access  to  a  shared 
communications  medium.  This  inherent  complexity  suggests  the  choice  of  a  MAC 
protocol  is  crucial  to  the  performance  of  a  WSN  [1].  There  are  many  options  for  MAC 
protocols,  and  the  best  choice  is  often  application  dependent.  Here  we  discuss  the  key 
features  and  benefits  of  the  IEEE  802.15.4  MAC  protocol  as  a  reference  framework  for 
MAC  operation. 

The  MAC  sub-layer,  which  appears  just  above  the  PHY  layer  of  the  OSI  model,  is 
responsible  for  managing  transmissions,  access  to  the  channel  and 
association/disassociation  to  the  network.  The  MAC  protocol  of  the  IEEE  802.15.4 
standard  supports  two  modes  of  transmission  [6]. 

•  Beacon-enabled:  In  this  mode,  periodic  beacons  are  transmitted  by  a 
designated  coordinating  node  (i.e.,  the  gateway)  in  order  to  maintain 
synchronization  and  exchange  network  information  between  the  nodes. 
These  beacons  constitute  the  first  slot  of  a  16  slot  superframe.  The 
exchange  of  data  takes  place  during  the  superframe.  A  superframe  is  a 
collection  of  timeslots  that  repeat  and  provide  a  format  for  scheduling 
transmissions.  The  structure  of  the  superframe  is  defined  by  the 
coordinating  node,  and  the  remaining  nodes  of  the  WSN  are  informed 
about  the  superframe  structure  through  beacons.  The  superframe  structure 

is  shown  in  Eigure  3.  It  can  be  seen  that  the  superframe  is  bounded  by  the 
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beacon  frames.  The  entire  superframe  may  be  occupied  by  the  contention 
access  period  (CAP).  During  the  CAP,  if  a  device  wishes  to  communicate, 
it  has  to  contend  with  other  devices  using  a  slotted  carrier-sense  multiple 
access  with  collision  avoidance  (CSMA/CA)  mechanism;  however,  if 
some  devices  (or  applications)  require  low  latency  communication,  then  a 
contention  free  period  (CFP)  is  introduced,  which  consists  of  guaranteed 
time  slots  (GTS)  [6].  The  inactive  period  allows  for  nodes  to  shut  off  for  a 
portion  of  the  communications  cycle  to  conserve  battery  power  and 
elongate  the  lifetime  of  the  WSN. 

•  Non  Beacon-enabled;  When  the  WSN  operates  in  this  mode,  there  are  no 
beacons  broadcast  by  the  coordinator  and  there  are  no  superframes. 
Access  to  the  channel  in  the  network  takes  place  using  an  unslotted 
CSMA/CA  mechanism  [6]. 


P(  me  It  o 
;  O  )S  s 
P(  ifi  0(  I 


Beacon  Frames 


Icjoi  it(  in  :ion| 


Inactive 

Period 


Figure  3.  IEEE  802.15.4  superframe  structure  from  [6]. 


Based  on  the  802.15.4  example,  it  can  be  seen  that  the  purpose  of  the  MAC  sub¬ 
layer  is  to  provide  coordinated  access  to  the  communications  channel  for  all  nodes  in  a 
WSN  in  an  effort  to  reduce  delay,  prevent  packet  collisions,  and  minimize  required  re¬ 
transmissions.  The  effectiveness  of  a  specific  MAC  implementation  plays  a  key  role  in 
the  power  consumption  of  the  nodes  and  is,  therefore,  a  critical  component  in  WSN 
design. 
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c. 


NETWORK  LAYER 


The  802.15.4  standard  does  not  address  the  staek  implementation  above  the  DLL 
and  instead  leaves  the  funetional  design  of  this  layer  to  the  network  engineer.  In  this 
section,  we  discuss  the  purpose  of  the  network  layer  and  the  considerations  that  need  to 
be  studied  when  designing  the  best  routing  algorithm.  The  network  layer  defines  how 
interconnected  networks  interact  with  one  another.  The  network  layer  is  the  first  layer  in 
the  OSl  model  that  is  concerned  with  sending  data  to  a  remote  network  in  contrast  to  the 
DLL,  which  functions  to  send  data  across  a  local  network  only  (i.e.,  communications  at 
the  DLL  are  point-to-point  rather  than  end-to-end).  Specific  functions  normally 
performed  by  the  network  layer  include; 

•  Logical  Addressing;  Every  node  that  communicates  within  a  WSN  must 
be  assigned  an  address.  A  familiar  example  of  addressing  is  the  Internet 
protocol  (IP)  address  for  the  network  layer  protocol.  Network  layer 
logical  addresses  are  independent  of  particular  hardware  and  must  be 
unique  across  an  entire  interconnected  network  [2]. 

•  Routing;  Moving  data  across  a  series  of  interconnected  networks  is  the 
defining  function  of  the  network  layer.  It  is  the  job  of  the  devices  and 
software  routines  that  function  at  the  network  layer  to  handle  incoming 
packets  from  various  sources,  determine  their  final  destination,  and  then 
determine  the  best  method  to  forward  them. 

•  Datagram  Encapsulation;  The  network  layer  encapsulates  messages 
from  higher  layers  by  placing  them  into  datagrams,  called  packets,  with  a 
network  layer  header. 

•  Fragmentation  and  Reassembly;  The  network  layer  must  send  messages 
down  to  the  DLL  for  transmission.  If  the  DLL  has  a  design  limit  on  the 
packet  length  that  differs  from  the  packet  length  used  at  the  network  layer, 
then  the  outgoing  packets  have  to  be  broken  into  fragments.  This  process, 
known  as  fragmentation,  must  be  un-done  when  packets  are  received  at 

the  destination  node.  This  is  known  as  reassembly. 
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•  Error  Handling  and  Diagnostics:  Special  protocols  are  used  at  the 
network  layer  to  allow  devices  that  are  logically  connected,  or  that  are 
routing  traffic,  to  exchange  information  about  the  status  of  hosts  on  the 
network  or  the  devices  themselves  [2], 

Network  layer  protocols  may  offer  either  connection-oriented  or  connectionless 
services  for  delivering  packets  across  the  network.  Connectionless  protocols  are  the  most 
common  at  the  network  layer.  In  many  protocol  suites,  the  network  layer  protocol  is 
connectionless,  and  connection-oriented  services  are  provided  by  the  upper  layers  [2]. 

WSN  specific  routing  requirements  and  energy  limitations  require  the 
development  of  routing  protocols  that  are  designed  for  the  efficient  and  effective 
operation  of  WSNs.  WSN  routing  protocols  can  be  divided  into  data-centric,  location 
based  and  hierarchical  types  of  routing  protocols.  In  data-centric  routing,  the  gateway 
sends  queries  to  certain  regions  of  the  network  and  waits  for  data  from  the  sensors 
located  in  the  selected  regions.  Sensor  Protocols  for  Information  via  Negotiation  was  the 
first  data-centric  protocol.  It  considers  data  negotiation  between  nodes  in  order  to 
eliminate  redundant  data  and  save  energy  [7].  Other  examples  of  notable  data-centric 
WSN  routing  protocols  are  gradient-based  routing  [8],  constrained  anisotropic  diffusion 
routing  [9],  and  active  query  forwarding  in  sensor  networks  [10].  Location-based  routing 
protocols  focus  on  delay  optimization  by  greedily  forwarding  packets  towards  sensors 
near  the  destination.  Location-based  protocols  include  minimum  energy  communication 
network  [11],  geographic  adaptive  fidelity  [12],  and  geographic  and  energy  aware  routing 
[13].  Hierarchical  WSN  routing  protocols  exploit  rigid  topology  constraints  to  provide 
energy-efficient  and  scalable  routing  [14].  Popular  hierarchical  routing  protocols  include 
the  low-energy  adaptive  clustering  hierarchy  [15]  and  power-efficient  gathering  in  sensor 
information  systems  [16]. 

D,  UPPER  LAYERS 

The  functions  of  the  layers  above  the  network  layer  in  the  OSI  model  include 
connection  establishment,  data  acknowledgement  and  retransmission,  data  multiplexing 
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and  de-multiplexing,  and  software  application  issues  [2],  The  functionality  of  the  upper 
layers  is  not  discussed  in  this  thesis.  Our  focus  is  on  the  PHY,  DLL,  and  network  layers 
of  the  WSN. 
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III.  INTERFERENCE  CONSIDERATIONS  IN  WIRELESS 
COMMUNICATIONS  SYSTEMS 


One  of  the  key  aspeets  that  must  be  eonsidered  when  modeling  interferenee  in 
wireless  eommunieations  is  how  interferenee  disturbs  the  reeeption  of  a  given  desired 
signal  [4],  Early  random-aeeess  models  assumed  that  if  any  two  transmitters  had 
temporally  overlapping  transmissions  then  the  transmissions  would  be  eollisions  and  all 
data  would  have  to  be  re-sent.  This  assumption  resulted  in  a  view  that  a  signal  would 
always  be  lost  in  the  presence  of  interference,  which  necessitated  the  use  of  medium 
access  protocols  that  provided  either  time  or  frequency  exclusivity.  Medium  access  of 
this  variety  provides  collision-free  operation  but  does  not  maximize  use  of  resources  by 
allowing  transmissions  from  nodes  with  sufficient  spatial  separation.  The  introduction  of 
the  radio  capture  concept  allows  for  the  successful  reception  of  a  desired  signal  in  an 
interference  environment  and,  thus,  provides  for  the  possibility  of  multiple  simultaneous 
transmissions.  These  transmissions  potentially  interfere  but  are  controlled  in  order  to 
allow  for  successful  reception  by  placing  design  limits  on  the  level  of  interference.  There 
are  many  radio  capture  models  available  in  the  literature,  which  can  be  broadly  grouped 
together  as  examples  of  protocol  or  physical  interference  models.  In  this  chapter  we 
develop  the  conceptual  basis  for  the  PrIM  and  the  PhIM  in  order  to  demonstrate  the 
superiority  of  the  PhIM  and  associated  trade-offs  of  using  the  more  accurate  model. 

A,  INTERFERENCE  MODEL  DEFINITIONS 

Consider  a  WSN  G(V,L)  consisting  of  n  wireless  nodes  F  =  {Vj, Vj,..., v„} , 

which  are  distributed  in  a  two-dimensional  planar  grid.  Each  node  is  wirelessly 

connected  to  its  neighboring  nodes  over  one  of  m  communications  links E  =  {/p/j,...,/^}  . 

In  the  development  of  the  PrIM  and  PhIM  in  [17],  arbitrary  and  random  network  settings 

are  provided  as  analysis  options.  In  the  arbitrary  network  setting,  the  location  of  nodes 

and  traffic  patterns  can  be  chosen  by  the  network  designer.  Since  any  traffic  pattern  and 

node  placement  can  be  used,  the  bounds  for  this  scenario  are  applicable  to  any  network; 

therefore,  the  arbitrary  network  setting  may  be  viewed  as  the  best  case  bounds  on 
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network  performance  [18].  In  the  random  setting  node  location,  traffic  patterns  and 
transmission  power  for  each  node  are  all  independently  selected  using  a  uniform  random 
distribution.  For  all  subsequent  analysis  in  this  thesis,  we  assume  arbitrary  node 
placement,  node  transmission  power,  and  traffic  pattern  settings  for  simplicity  of 
analysis. 

B,  PROTOCOL  INTERFERENCE  MODELS 

Under  the  protocol  model,  the  impact  of  interference  from  neighboring  nodes  is 
binary  and  is  solely  determined  by  whether  or  not  the  node  falls  within  the  interference 
range  of  non- intended  transmitters  [19].  The  conceptual  basis  for  the  protocol 
interference  model  was  initially  introduced  in  [17]  in  an  effort  to  analyze  the  capacity  of 
multi-hop  wireless  networks.  The  protocol  interference  model  has  been  the  focus  of 
rigorous  study  and  application  since  being  introduced  and  is  the  basis  for  numerous  MAC 
protocol  implementations.  The  basis  of  the  protocol  interference  model  is  the 
vulnerability  capture  model,  which  states  that  the  receiving  node  that  is  closest  to  the 
transmitting  node  in  a  given  wireless  network  has  the  greatest  received  power. 

The  PrIM  is  an  implicit  interference  model  that  describes  the  effects  of 
interference  based  on  the  communications  channel  conditions  [4].  Given  the  WSN 
G(V,  L) ,  the  PrIM  asserts  that  a  transmission  from  y  to  v j ,  located  at  Euclidean 

distance  d  (i,  j) ,  over  some  link  eL,  is  successful  if  for  every  other  node  that  is 
simultaneously  transmitting  is  located  d (k,  j)  away.  This  is  stated  as 

d{k,j)>{l  +  A)d{i,j),A>0  (3.1) 

where  A  is  a  guard  parameter  to  ensure  that  all  conflicting  transmitting  nodes  are 
sufficiently  far  away  from  the  receiver  to  allow  successful  decoding  of  the  intended 
signal  [18].  In  order  to  relate  the  A  parameter  to  a  physically  significant  quantity,  the 
PrIM  introduces  the  concept  of  the  maximum  transmission  range  and  the  maximum 
interference  range  where  and  represent  the  maximum  transmission  and 

maximum  interference  ranges  of  each  node,  respectively.  The  transmission  range 
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represents  the  maximum  distanee  up  to  whieh  a  packet  can  be  received,  while  the 
interference  range  represents  the  maximum  distance  up  to  which  simultaneous 
transmissions  interfere.  In  the  literature,  the  interference  range  is  usually  chosen  to  be 
twice  as  large  as  the  transmission  range,  which  has  been  shown  to  not  necessarily  be  a 
practical  assumption  [5].  In  reality,  the  actual  values  of  the  transmission  and  interference 
ranges  depend  on  the  transmission  power  used  by  the  nodes  and  are  a  function  of  the  SIR 
threshold  [5].  Optimized  and  more  realistic  values  for  can  be  obtained  using  a 
reality  check  mechanism  as  shown  in  [19],  which  helps  to  obtain  a  more  realistic 
agreement  between  the  protocol  model  behavior  and  real-world  physical  interference 
behavior. 

Within  the  construct  of  the  PrIM,  there  are  two  different  types  of  interference  that 
have  been  studied  in  the  literature,  namely  primary  interference  and  secondary 
interference.  Primary  interference  occurs  when  a  node  transmits  and  receives  packets  at 
the  same  time.  Secondary  interference  occurs  when  a  node  receives  two  or  more  separate 
transmissions  [20]. 

The  inadequacies  of  the  PrIM  become  apparent  when  analyzing  the  secondary 
interference  problem.  The  assumption  that  concurrent  transmissions  do  not  interfere  as 
long  as  the  nodes  that  are  exchanging  information  are  outside  of  one  another’s 
interference  radius  does  not  accurately  account  for  aggregate  interference.  According  to 
the  PrIM,  the  aggregate  interference  stems  from  the  fact  that  multiple  simultaneous 
transmissions  with  sufficient  spatial  separation  can  combine  to  constructively  interfere 
and  prevent  successful  reception. 

C.  PHYSICAL  INTERFERENCE  MODELS 

The  PhIM  is  based  on  the  power  capture  model  that  assumes  a  threshold  based 
channel  model.  In  other  words,  the  transmission  data  rate  is  a  function  of  the  signal-to- 
interference  ratio  (SIR).  If  the  SIR  is  sufficiently  large,  the  data  rate  is  guaranteed  and 
constant.  Conversely,  if  the  SIR  is  below  a  pre-defined  threshold,  then  the  transmission 
data  rate  is  zero  [4].  The  SIR  is  calculated  from 
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(3.2) 


SIR  = 
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where  P  is  the  SIR  threshold  above  which  data  transmissions  are  successful  from  to 
Vj  in  the  presence  of  all  concurrent  secondary  interference  events,  V  is  the  subset  of 
nodes  in  the  WSN  that  are  assigned  to  transmit  simultaneously  and  Py(/)  is  the  received 
power  at  node  j  due  to  node  i.  When  calculating  the  SIR,  P.{i) is  calculated  from 


P.{i)  =  d{i,jr\a>2  (3.3) 

where  d  (i,  j)  is  the  Euclidian  distance  between  nodes  and  a  is  the  path  loss  exponent, 
which  is  a  function  of  the  propagation  loss  model  used  to  account  for  free-space  power 
loss  during  radio  transmission. 

As  can  be  seen  from  the  above  equations,  the  PhIM  does  not  suffer  from  the  same 
shortcomings  as  the  PrIM.  Namely,  the  error  introduced  by  neglecting  aggregate 
interference  effects  and  the  unrealistic  binary  channel  condition  utilized  in  the  PrIM  are 
not  concerns  when  using  the  PhIM.  These  advantages  in  accuracy  are  countered  by 
significantly  increased  complexity  in  simulation  analysis.  Additionally,  the  complexity  of 
PhIM  based  MAC  protocol  implementations  is  greatly  increased  as  compared  to  available 
PrIM  solutions. 


D,  GRAPH  BASED  INTERFERENCE  MODELS 

Graph  theory  is  an  important  analytical  tool  in  the  field  of  wireless  network 
design  and  modeling.  Graph  theory  provides  the  ability  to  describe  a  network  with  a  set 
of  vertices  and  edges  that  connect  the  vertices.  A  graph  G  is  generally  defined  as  a 
function  of  its  vertices  and  edges  (i.e.,  G(V,E)).  In  a  WSN,  the  sensor  nodes  are  the 
vertices,  and  the  wireless  links  between  nodes  are  the  edges.  An  example  connectivity 
graph  is  shown  in  Figure  4  with  a  set  of  vertices  V  =  {A,B,C,D,E]  connected  by  a  set  of 

edges  E  =  [AB,AE,BC,BD,BE,CD,DE]  .  The  graph  G  is  said  to  be  connected  if  there 
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is  a  path  from  any  vertex  to  any  other  vertex.  Graphs  are  either  undirected,  in  which 
edges  are  represented  by  a  line,  or  directed  (also  called  a  digraph),  in  which  an  edge  is 
represented  by  an  arrow  drawn  from  the  tail  of  a  vertex  to  the  head  of  a  vertex  [4] . 

G  =  (V,E) 


A 


Figure  4.  Example  of  an  undirected  graph  from  [8]. 

A  basic  aspect  of  an  ad-hoc  wireless  network  that  can  be  modeled  by  means  of  a 
graph  is  the  radio  connection  between  any  two  nodes.  A  typical  rule  to  determine  if  a 
given  node  v.  can  communicate  with  another  node  Vj  is  that  the  SIR  of  the  received 

signal  at  v-  must  be  equal  to  or  larger  than  a  given  threshold  (3  [8].  If  this  connectivity 
criterion  is  satisfied,  then  an  edge  is  assigned  from  v.  to  Vj .  Thus,  a  connectivity  graph 

representing  a  wireless  network,  such  as  a  WSN,  can  be  created  with  knowledge  of  node 
location  and  node  transmit  power.  The  network  connectivity  described  by  a  connectivity 
graph  can  be  misleading,  and  it  should  be  noted  that  this  connectivity  does  not  account 
for  interference  due  to  concurrent  transmissions.  The  interference  encountered  in 
networks  that  utilize  a  packet  forwarding  scheme,  allowing  simultaneous  transmissions, 
is  usually  accounted  for  by  constructing  a  structure  known  as  a  conflict  graph  [5].  In  a 
conflict  graph,  denoted  as  G, ,  a  vertex  is  introduced  for  each  link  in  the  original 

network  G{V,E)  .  An  edge  in  (F, ,  )  connects  two  vertices  in  the  conflict  graph  if 

these  two  links  interfere  in  the  corresponding  connectivity  graph  (as  in  the  case  of 
simultaneous  transmissions). 
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The  appropriate  parameter  of  a  wireless  link  to  eharacterize  the  interference 
relationship  between  two  links  is  the  SIR  tolerated  at  the  receiver.  At  this  point,  we 
distinguish  two  approaches  for  assigning  values  to  the  links  in  G;{Vj,Ej):  (1)  in  its 
aggregate  form,  as  in  the  PhIM,  or  (2)  by  considering  interference  from  each  interferer 
individually,  as  in  the  PrIM  [4],  [5].  As  stated  in  the  previous  discussion,  the  PrIM  has  a 
and  associated  with  each  node  that  can  be  determined  by  using  a  free  space 

path-loss  model  and  the  node  transmit  power.  In  the  case  of  a  PrIM  based  conflict  graph, 
an  edge  is  added  to  the  graph  in  the  event  of  secondary  interference  by  a  simultaneously 
transmitting  node.  Note  that  the  interference  between  two  links,  according  to  this  rule,  is 
not  a  reciprocal  relationship,  and  interference  graphs  are  in  general  directed  [4]. 

Recalling  that  the  PrIM  is  a  binary  communications  model,  for  every  instance  that 
a  link  exists  in  the  conflict  graph,  we  see  that  there  is  a  need  for  exclusivity  in  channel 
access  for  the  conflicting  nodes.  Thus,  the  number  of  nodes  permitted  to  transmit  at  the 
same  time  (i.e.,  the  amount  of  link  reuse)  is  a  direct  output  of  the  conflict  graph.  The 
amount  of  link  reuse  allowed  is,  therefore,  a  function  of  the  loss  propagation  model, 
transmission  power,  and  interference  range  model  used.  In  contrast  to  the  PrIM,  the 
PhIM  is  an  aggregate  interference  model;  thus,  the  conflict  graph  link  weights  consist  of 
aggregate  SIR  values  for  all  simultaneously  transmitting  nodes.  The  PhIM  takes  into 
account  all  active  links  regardless  of  the  distance  between  the  links,  which  results  in  a 
distinctly  different  conflict  graph  as  compared  to  the  PrIM  for  a  given  network  topology. 
The  PhIM  conflict  graph  presents  a  more  accurate  representation  of  the  interference  in 
the  physical  environment  as  compared  to  the  PrIM.  Unfortunately,  in  general  it  is 
difficult  to  obtain  accurate  and  continuously  updated  SIR  information;  additionally,  there 
is  an  overhead  incurred  to  share  that  SIR  information  for  dynamically  assigning  conflict 
graph  link  weights.  This  extra  computation  and  the  associated  power  and  bandwidth 
consumption  results  in  the  choice  of  a  PrIM  approach  to  interference  modeling  for  many 
models.  In  this  thesis  we  propose  a  simplified  approach  to  assigning  SIR  values  to  a 
PhIM  conflict  graph  in  an  effort  to  bridge  the  gap  between  the  simplicity  of  the  PrIM  and 
the  accuracy  of  the  PhIM. 
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IV.  MEDIUM  ACCESS  CONTROL 


In  this  chapter  we  shift  our  focus  from  the  accuracy  and  computational  ease  with 
which  interference  can  be  modeled  in  a  WSN  to  how  we  use  our  interference  analysis 
results  to  best  share  the  communications  resources  amongst  the  WSN  nodes.  The  MAC 
protocol  is  primarily  responsible  for  orderly  and  efficient  resource  management  in  a 
wireless  networking  system.  The  medium  access  meehanism  can  be  implemented  in 
using  either  a  distributed  coordination  function  or  a  centralized  coordination  function. 
Generally  speaking,  a  distributed  coordination  function  consists  of  a  contention  based 
medium  access  method  where  nodes  compete  with  their  neighbors  for  transmission  time 
and  channel  resources.  Central  coordination  functionality,  on  the  other  hand,  involves 
the  execution  of  a  carefully  synchronized  transmission  schedule,  which  is  determined  by 
a  central  entity.  An  example  of  a  coordinated  MAC  protocol  is  the  time-division  multiple 
access  (TDMA)  protocol,  which  allocates  network  resources,  node  by  node,  using  time 
slots.  There  are  inherent  complexities  associated  with  using  TDMA  including  timing 
synchronization,  message  exchange  for  slot  assignment,  and  slot  schedule  optimization  to 
account  for  various  network  dynamics.  In  this  chapter  we  present  a  detailed  explanation 
of  the  TDMA  protocol  and  a  discussion  of  the  advantages  of  using  TDMA  for  a  WSN  as 
opposed  to  a  contention  based  MAC  protocol. 

An  upgrade  to  TDMA  is  the  STDMA  protocol,  which  takes  advantage  of  the 
spatial  separation  of  nodes  in  a  WSN  and  allows  more  than  one  node  to  transmit  in  the 
same  time  slot  (slot  reuse).  STDMA  is  superior  to  TDMA  and  can  decrease  delay  and 
increase  network  capacity.  In  order  for  an  STDMA  protocol  to  operate  properly,  there 
must  be  a  thorough  understanding  of  the  interference  environment  created  by  allowing 
the  time  slot  reuse.  In  this  chapter  we  explain  the  operation  of  the  TDMA  and  STDMA 
protocols,  compare  and  contrast  the  two  protocols,  demonstrate  the  superiority  of 
STDMA  to  TDMA,  and  discuss  the  effects  of  accurate  interference  modeling  on  the 
performance  of  an  STDMA  protocol  implementation. 
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A,  TIME-DIVISION  MULTIPLE  ACCESS 

In  a  TDMA  access  scheme,  a  conflict-free  access  schedule  is  aehieved  by 
assigning  a  time  slot  to  eaeh  node  and  only  one  node  can  transmit  per  time  slot.  A  frame 
is  formed  with  the  concatenation  of  all  time  slots  and  the  required  non-data  overhead 
timeslots.  In  a  traditional  TDMA  protoeol,  each  node  is  assigned  one  time  slot  per  frame, 
and  the  slots  are  assigned  by  a  controlling  entity,  such  as  the  gateway  in  a  WSN.  The 
time  slot  schedule  ean  be  predefined  to  aid  in  simplicity  of  network  design  or  it  ean  be 
periodically  updated  to  refieet  changes  in  network  traffie  patterns  to  increase 
performanee.  In  this  thesis,  we  use  the  simplified  TDMA  protoeol  with  a  pre-defined  slot 
assignment  sehedule  as  our  reference  TDMA  model. 

In  a  TDMA  seheme  the  time  axis  is  divided  into  M  time  slots  (assuming  a  WSN 
consisting  of  M  nodes)  that  are  pre-assigned  to  the  different  WSN  nodes.  Eaeh  node  is 
permitted  to  transmit  using  the  entire  channel  bandwidth  during  its  assigned  time  slot. 
The  slot  assignments  follow  a  predetermined  pattern  that  repeats  itself  periodically,  and 
each  such  period  is  called  a  frame.  In  the  most  basic  TDMA  scheme,  every  user  has 
exaetly  one  slot  in  every  frame  as  indieated  by  the  numerieal  indices  in  Figure  5. 

In  our  analysis  of  TDMA,  we  are  interested  primarily  in  the  throughput  and  delay 
eharaeteristies,  which  depend  on  frame  length,  slot  assignment,  and  slot  assignment 
ordering.  We  consider  the  throughput  to  be  the  average  aggregate  data  that  ean  be 
transported  through  the  channel  per  unit  time.  We  quantify  the  throughput  as  the  fraction 
of  time  during  whieh  node  data  is  being  suceessfully  transferred.  The  delay  is  the  time 
from  message  generation  until  message  delivery  to  the  next  hop.  We  infer  that  a  lower 
per-hop  delay  indicates  a  lower  overall  system  wide  delay  when  sending  information 
from  an  arbitrary  node  to  the  gateway. 
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Figure  5.  Example  of  a  TDMA  frame  structure  from  [21]. 
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To  analyze  the  performanee  of  a  TDMA  seheme,  eonsider  a  system  eomposed  of 
M  nodes  with  eaeh  node  transmitting  equally  long  paekets  of  P  bits  eaeh.  If  the  total 
rate  of  transmission  is  R  bits/sec,  then  the  paeket  transmission  time  is 


T  = 


P 

—sec . 
R 


(4.1) 


This  is  taken  to  be  the  slot  size.  The  duration  of  the  entire  cycle  is,  therefore,  equal  to 
where 


T^=MTsqc.  (4.2) 

Assuming  that  the  packet  generation  processes  of  new  packets  by  each  WSN 
node  are  independent,  it  follows  that  the  queuing  behaviors  of  each  nodes’  queue  are 
independent.  This  assumption  that  the  queues  are  independent  is  reasonable  because  each 
node  transmits  a  packet  every  7)  seconds,  independently  of  any  event  in  any  of  the  other 

queues  of  other  users.  Consequently,  in  the  following  analysis,  we  concentrate  on  the 
characteristics  of  one  node.  Without  loss  of  generality,  we  assume  that  the  node 
transmits  a  packet  if  there  is  a  packet  in  the  node’s  queue  during  its  assigned  slot  of  every 
frame. 

Consider  a  typical  packet  generated  by  a  node.  The  delay  suffered  by  this  packet 
has  three  components:  (1)  the  time  between  its  generation  and  the  end  of  the  current 
frame,  (2)  the  queuing  time  to  allow  all  the  packets  already  queued  to  be  transmitted  and, 
(3)  the  packet  transmission  time  itself.  Of  these  components,  the  first  and  the  third  are 
readily  known  since  all  frames  are  of  equal  length,  and  the  average  time  between  the 
packet  generation  time  and  the  end  of  the  current  frame  is  0.57^ .  Here  we  assume  that  the 

average  packet  is  generated  halfway  through  an  arbitrary  frame  period  (i.e.,  0.5  of  Tj- 
By  restricting  our  packet  generation  to  Poisson  arrivals  and  using  the  deterministic 
servicing  of  TDMA,  our  queuing  time  for  each  node  can  be  described  by  an  M/D/1  queue 


w  = — ^ — 

•  2(l-p) 


X 


(4.3) 
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where  x  -T^  is  the  deterministic  service  time  and  p^  =  is  the  fraction  of  the  total 
network  load  arriving  at  node  i  [22],  The  packet  generation  at  the  node  is  A.  packets  per 
second.  Note  that  the  queue  notation  M/D/1  denotes  a  single  server  queue  with  Poisson 
arrivals  and  a  deterministic  service  time.  Substituting  for  X  and  recalling  that 

=  MT  ,  we  obtain  the  expected  queuing  time 


^ — 
^  2(1 -p) 


MT  sec . 


(4.4) 


Combining  Eqs.  (4.2)  and  (4.3),  for  an  average  packet  generated  at  time  0.57)  within  a 
frame  and  with  a  transmission  time  of  T  sec  ,  we  obtain  the  total  packet  delay  [22] 


1  ‘  ^  2  2(1-/?) 


MT+T=T 


1  + 


M 


2(l-p) 


sec . 


(4.5) 


From  [22]  the  relationship  between  throughput  S  and  utilization/?  for  a  M/D/1 
queue  is  shown  to  be  S  =  p  .  Substituting  S  for  p  in  Eq.  4.5,  we  get 


D  ,  M 

—  —  1  ■! - 7 - r 

T  2(1-5) 


(4.6) 


Equation  4.6  is  the  TDMA  delay  throughput  characteristic  where  M  is  the  number  of 
nodes  in  the  network,  S  is  the  network  throughput,  and  D  is  the  delay  in  the  network.  A 
family  of  graphs  of  the  expected  delay  versus  throughput  for  various  network  sizes  (i.e., 
different  values  of  M  )  is  shown  in  Figure  6. 

Though  TDMA  provides  for  simple  implementation  as  well  as  closed  form 
equations  to  be  used  in  design  and  analysis,  it  does  not  provide  for  optimal  utilization  of 
resources  on  a  per  node  basis.  This  drawback  is  particularly  harmful  in  the  context  of  an 
operating  WSN  where  battery  power  is  limited  and  care  must  be  taken  to  maximize 
resource  utilization.  In  the  following  section,  we  discuss  how  resource  utilization  can  be 
improved  by  allowing  simultaneous  transmissions  for  multiple  nodes  while  mitigating 
inter-node  interference  [23]. 
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Figure  6.  TDMA  throughput  versus  expeeted  delay  curve  for  various  network  sizes 

from  [22], 

B,  TIME-DIVISION  MULTIPLE  ACCESS  WITH  SPATIAL  LINK  REUSE 

Access  schemes  based  on  the  controlled  simultaneous  transmissions  for  nodes 
with  sufficient  spatial  separations  are  known  as  STDMA  protocols.  STDMA  operates  by 
allowing  multiple  nodes  to  transmit  at  once  as  long  as  the  nodes  are  sufficiently  separated 
in  space.  In  this  case  a  given  node  may  have  the  opportunity  to  transmit  more  than  once 
per  frame,  thereby  decreasing  the  average  delay  at  each  node.  In  order  for  this  delay 
reduction  to  be  useful,  it  is  important  that  there  is  a  thorough  understanding  of  the 
interference  environment  created  by  using  STDMA.  An  adequate  model  must  be  used  to 
represent  network  interference  in  order  to  predict  network  performance  under  STDMA. 

Similar  to  TDMA,  STDMA  is  a  synchronized  protocol  that  uses  a  centralized 
coordination  approach.  STDMA  requires  an  algorithm  for  generating  the  transmission 
schedule  [23].  Transmission  scheduling  using  STDMA  is  a  complex  problem  that 
requires  optimization  in  order  to  maximize  resource  utilization.  Ideal  characteristics  of 
an  STDMA  scheduling  algorithm  are  as  follows:  (1)  the  algorithm  must  be  able  to  adapt 
to  changes  in  network  topology  and  (2)  the  algorithm  must  be  simple  enough  to  allow  for 
fast,  on-line  computation.  In  addition  to  computing  transmissions  schedules  that  allow 
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for  operation  in  the  presence  of  interference,  the  STDMA  schedule  may  be  further 
optimized  to  dynamically  re-order  time  slots  in  reaction  to  changing  traffic  patterns  [24], 
In  practice  there  is  no  optimal  STDMA  algorithm  due  to  the  computational  complexity  of 
the  STDMA  implementation;  therefore,  the  common  approach  is  to  develop  fast  heuristic 
algorithms  to  solve  the  STDMA  problem  [23]. 

In  the  remainder  of  this  section,  our  goal  is  to  show  the  derivation  of  the 
maximum  throughput  for  an  STDMA  protocol  assuming  knowledge  of  the  interference 
environment  and  that  there  is  no  time  constraint  for  scheduling  computation.  It  is 
important  to  note  that  when  scheduling  a  transmission  in  an  STDMA  protocol,  the  time 
slots  can  be  assigned  to  either  nodes  or  links.  In  a  node  oriented  assignment,  a  node  that 
is  scheduled  to  transmit  may  transmit  to  any  other  node,  which  allows  for  broadcast  and 
multicast  transmissions.  For  link-based  transmission  schedules,  each  node  is  only 
permitted  to  transmit  over  a  specified  point-to-point  link  to  another  node.  Clearly,  the 
broadcast  and  multicast  capability  afforded  by  node-based  scheduling  results  in  a  higher 
number  of  available  active  links  per  node  as  compared  to  a  link-based  schedule.  This 
characteristic  of  node-based  STDMA  serves  as  motivation  for  us  to  use  a  node-based 
STDMA  schedule  in  our  presentation  of  the  comparison  of  STDMA  to  TDMA 
throughput. 

For  the  purposes  of  the  following  discussion,  our  network  model  consists  of  a 
connectivity  graph  G(V,E)  and  utilizes  the  PrIM  for  describing  connectivity  and 
interference.  We  assume  that  there  are  directional  links  between  nodes,  there  are  no 
primary  interference  conflicts  allowed,  and  each  link  is  error  free  as  long  as  the  SIR  >  P  . 
We  assume  all  nodes  use  the  same  transmit  power,  and  we  use  frame  and  timing  notation 
consistent  with  our  earlier  discussions  on  TDMA  and  SIR. 

Recall  that  2.  is  the  average  traffic  load  arriving  at  node  i  and  let  /?.  denote  the 

number  of  slots  in  which  node  i  is  scheduled  using  a  node-oriented  schedule  [23].  The 
throughput  is  the  maximum  A.  for  which  the  delay  remains  finite,  which  requires  the 
following  inequality  to  hold; 
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(4.7) 


The  throughput  is  the  largest  A-  that  satisfies  Eq.  4.7.  This  is  represented  as 


^  ieN  p  /sec 


(4.8) 


where  p.  =  — 
'  A 


[23]. 


Having  provided  this  formal  definition  of  throughput,  our  goal  is 


now  to  maximize  the  throughput  for  a  given  node  assignment  schedule.  We  proceed  by 

introducing  the  concept  of  a  transmission  group,  which  is  a  group  of  nodes  that  can  share 

L  s 

a  time  slot.  Let  ^  denote  the  sets  of  transmission  groups  and  "  serve  as  a  binary 

indication  parameter  defined  as  [23] 


[  1,  node/e/ 
[O,  otherwise 


(4.9) 


Let  X;  be  the  number  of  time  slots  assigned  to  transmission  group  /.  The 
maximization  of  the  throughput  z  is 


I, 


- ,V/eZ 


PT. 

subject  to  the  following  constraints  [23]; 


'N 


X,  >  0,  X;  e  □  ,  V/  e 


(4.10) 


(4.11) 

(4.12) 


where  □  is  the  set  of  integers. 

There  exists  a  sizeable  field  of  study  dedicated  to  providing  optimal  and  fast 
solutions  to  this  optimization  problem.  The  solution  of  this  optimization  requires  the  use 
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of  linear  programming  numerical  methods.  The  results  of  the  optimization  are  detailed  in 
[23]  and  show  that  the  throughput  for  STDMA  can  be  as  much  as  35  times  the  throughput 
in  a  traditional  TDMA  system. 
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V.  INTERFERENCE  AWARE  LEAST-COST  ROUTING 


In  this  chapter  we  present  our  proposed  method  of  improving  upon  the 
performanee  of  traditional  TDMA  by  implementing  a  physieal  interferenee  based 
network  routing  algorithm.  Our  method  consists  of  a  SIR  threshold  based  PHY  layer 
reception  model,  STDMA  DLL,  and  shortest  path  routing  at  the  network  layer.  The 
discussions  on  interference  modeling,  graph  theory  and  STDMA  in  the  preeeding 
ehapters  provide  the  background  in  explaining  our  proposed  WSN  eommunications  stack. 

There  are  various  algorithms  for  finding  the  shortest  path  if  the  edges  in  a  network 
are  restrieted  to  non-negative  values.  The  most  popular  shortest  path  algorithm  in  use  is 
Dijkstra’s  algorithm  [25].  The  operation  of  Dijkstra’s  algorithm  in  the  eontext  of  using 
SIR  values  as  the  link-eost  metric  associated  with  transmission  groups  in  an  STDMA 
MAC  protocol  is  explained  in  this  chapter.  We  propose  to  use  the  PhIM  to  develop  a 
network  conflict  graph  in  order  to  quantify  the  interferenee  encountered  at  a  given 
network  node  during  data  exehange  over  a  wireless  link.  Once  these  interference  values 
are  obtained  and  link  weights  assigned,  we  use  Dijkstra’s  algorithm  to  find  the  shortest 
path  across  a  WSN  using  the  interferenee  weights  as  the  routing  metrie.  In  Chapter  VI, 
we  simulate  this  routing  method  in  conjunetion  with  a  STDMA  MAC  protoeol  and 
demonstrate  that  using  interferenee-aware  least  cost  routing  is  an  effective  means  of 
redueing  network  delay  when  compared  to  interferenee  aware  routing  using  the 
traditional  TDMA  MAC  protocol. 

A,  DIJKSTRA’S  ALGORITHM 

The  problem  of  finding  shortest  paths  plays  a  central  role  in  the  design  and 
analysis  of  networks.  Most  routing  problems  can  be  solved  as  shortest  path  problems 
once  an  appropriate  cost  is  assigned  to  each  link.  For  example,  eost  metries  can  reflect 
available  bandwidth,  delay  or  interferenee  [25].  In  a  packet  switehed  communieations 
network  sueh  as  a  WSN,  the  funetion  of  routing  data  across  a  network  using  the  shortest 
path  is  aecomplished  at  the  network  layer. 
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Dijkstra’s  algorithm  finds  the  least-cost  distance  from  any  node  to  any  other  node 
in  a  directed  graph  with  non-negative  weights.  Dijkstra’s  algorithm  has  been  shown  to 

work  in  time  for  a  network  with  n  nodes  with  reduced  complexity  in  the  instance 

of  sparse  connectivity  (i.e.,  non-dense  networks).  The  inputs  to  the  algorithm  are  the 
graph  G(k,£') ,  weights  c :  £'(G) and  a  starting  node 5  e  F(G)  where  E[G)  and 

V (G)  denote  the  node  and  link  weights  of  G(V,E)  ,  respectively.  The  output  consists 
of  the  shortest  path  from  s  to  all  v  e  F  (G) ,  and  the  associated  path  lengths  are  denoted  as 
the  sets  /(v)  and (v) ,  respectively. 

The  steps  of  the  algorithm  are  as  follows  [26]: 

1)  Set  /(5)  □  0 ,  set  / (v)  □  co  for  all  v  e  F (G)  \  {5} ,  set  R:=0 

2)  Find  a  vertex  j  v  e  F(G)\i? : /(v)  =  min[/(w)] 

[  wer(G)\R 

3)  Set 

4)  For  all  |we  F(G)\i? :  (v,w)  e  F'(G)|  do: 

If  /(w)  >/(v)  +  c((v,w))  then 

set  /(w)  := /(v)-i-c((v, w))  and  p(w):=v. 

5)  If R  F(GJ  then  go  to  2) 

Step  one  in  the  algorithm  sets  node  s  as  the  starting  node  and  assigns  a  cost  of 
zero  to  node  s  while  assigning  an  infinite  initial  cost  to  all  other  nodes.  Additionally,  the 
set  of  shortest  paths  is  initialized  as  empty.  In  step  two,  the  shortest  path  to  all  nodes 
from  the  current  node  is  found  and  this  path  is  added  to  the  set  R  in  step  three.  In  step 
four  all  paths  are  compared  to  the  newest  visited  node  and  the  current  path  is  verified  as 
the  shortest  path.  If  this  current  path  is  not  the  shortest  path,  then  the  shortest  path  is 
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calculated  based  on  the  current  cost  value.  This  shortest  path  overwrites  the  current  path. 
Verification  that  all  nodes  have  been  visited  and  termination  of  the  algorithm  once  all 
nodes  have  been  visited  takes  plaee  in  step  5. 

B,  INTERFERENCE  AWARE  LEAST-COST  ROUTING  ALGORITHM 

Dijkstra’s  algorithm  has  been  proven  to  find  the  shortest  path  among  all  nodes  in 
a  network  and  is,  therefore,  a  logical  choice  for  finding  the  shortest  path  through  a  WSN. 
In  order  for  us  to  integrate  Dijkstra’s  algorithm  with  our  overall  task  of  finding  least 
interfering  paths,  we  must  first  construct  the  eonneetivity  and  cost  graphs  for  our 
STDMA  WSN  scenario.  A  depiction  of  the  physical  layout  of  a  21-node  WSN  is  shown 
in  Figure  7.  Each  node  is  plaeed  in  a  Cartesian  coordinate  system  in  10-meter  inerements. 
We  ereated  the  network  shown  in  Figure  7  using  MATLAB  [27].  The  gateway  for  this 
WSN  is  designated  as  node  21,  and  all  other  nodes  are  considered  sensing  nodes  that  send 
information  to  the  gateway.  The  set  of  nodes  corresponding  to  this  speeific  WSN  model 
is  E(G)  =  {1,2,3,...,21}. 


Network  topology 
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Figure  7.  21 -Node  WSN  topology  generated  using  MATLAB,  where  node  21  is  the 

designated  gateway. 
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Given  this  topology,  we  construet  the  conneetivity  graph  starting  with  the  allowed 
conneetions  shown  in  Figure  8.  Allowed  connections  are  defined  as  connections  that 
forward  information  toward  the  gateway. 
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Figure  8.  Connectivity  graph  for  the  21-node  WSN. 

For  our  analysis,  the  term  “links”  refers  to  the  directed  edges  that  belong  to  the  set 
of  allowed  connections.  The  set  of  allowed  connections  is  denoted  as  L{G) .  The  links  for 

this  connection  scenario  areT(G)  =  A-6’A-7’-"’A7-2i’A8-2i’ A9-21}  •  Now  that  we  have 

obtained  V (G)  and  T(G) ,  we  use  these  two  sets  to  compose  graph  G(V,L) . 
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The  next  step  is  to  take  our  knowledge  of  the  topology  and  compute  the  power 
received  at  each  node  assuming  10  mW  transmission  from  a  node  that  corresponds  to  the 
maximum  transmit  power  of  a  Mica2  mote,  which  is  a  common  commercial  sensor  node 
[28].  We  compute  this  received  power  from 

P,{j)  =  P^(d(i.j)f  (5.1) 

where  Pi{j)  represents  the  power  received  at  node  j  due  to  a  transmission  by  node  1, 

=\0  mW  ,  d  (l,y)  is  the  distance  from  node  1  to  node  j,  and  a  =  2  is  the  path  loss 
factor.  This  calculation  must  be  repeated  for  all  nodes  where  yV  1  and  repeated  again 
for  nodes  1  through  20.  Since  node  21  is  the  gateway,  there  are  no  transmissions  from 
this  node;  thus,  there  is  no  requirement  to  compute  the  received  power  at  the  other 
network  nodes  due  to  transmissions  from  node  21 .  These  received  power  values  are  used 
to  compute  the  aggregate  interference  in  our  PhIM  conflict  graph  assignments. 

Given  our  goal  of  maximizing  spatial  reuse  and  controlling  interference,  we  do 
not  allow  one-hop  neighbors  to  simultaneously  transmit  because  of  the  creation  of  excess 
interference.  We  allow  for  two-hop  neighbors  to  transmit,  which  results  in  the  four  color- 
coded  transmission  groups  shown  in  Figure  9.  The  colors  denote  nodes  that  are  permitted 
to  simultaneously  transmit  in  our  STDMA  schedule.  For  example,  all  blue  nodes  may 
transmit  simultaneously  in  the  same  time  slot.  This  holds  for  all  the  colored  node  groups. 
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Network  topology 


Figure  9.  Transmission  groups  in  the  21-node  WSN  to  indieate  simultaneous 
transmissions  in  the  STDMA  transmission  sehedule. 

The  next  step  is  to  eompute  the  interference  experienced  during  each  reception 
due  to  simultaneous  transmissions  throughout  the  transmission  group.  For  example,  what 
interference  does  node  2  experience  when  listening  to  node  1  with  the  interference  due  to 
transmissions  from  nodes  3,  5,  11,  13,  and  15?  Note  that  nodes  1,  3,  5,  11,  13,  and  15  all 
belong  to  the  same  transmission  group.  For  this  calculation  we  assume  that  all  nodes  in 
the  transmission  group  have  data  to  send,  which  provides  the  maximum  possible 
interference.  The  calculations  for  the  SIR  at  each  node  are  calculated  from 

keTGi 

where  TG-  denotes  the  transmission  group  of  node  j  .  This  was  also  given  in  Eq.  (3.2). 

This  states  that  the  interference  received  at  node  i  during  reception  from  transmitting 
node  j  (i.e.,  node  2  from  node  1)  is  the  aggregate  of  the  transmitted  power  from  all  other 
nodes  in  the  transmission  group  of  node  j  .  Equation  (5.2)  maintains  the  assumption  that 
all  nodes  in  the  transmission  group  send  data  simultaneously  and  that  the  interference  is 
at  its  maximum  possible  value  (upper  bound  on  experienced  interference).  This 
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calculation  is  repeated  for  every  combination  of  receiving  nodes  and  applieable 
transmission  groups.  This  enables  us  to  assign  an  interference  metric  (i.e.,  the  calculated 
SIR  value)  to  each  edge  in  our  interference  graph. 


To  ereate  the  confliet  graph,  we  treat  the  SIR  values  from  the  previous  seetion  as 
link  weights  such  that  higher  SIR  values  equate  to  a  lower  cost  for  transmission. 
Specifically,  we  assign  to  eaeh  link  in  T(G)  the  value  of  the  inverse  SIR  for  the 

receiving  node  of  the  link.  For  example,  the  link  weight  for  /j  j  is  assigned  related  to  the 
SIR  value  corresponding  to 


n-2 


[5ffi,(l)] 


f.(l) 

E  A(u 


(5.3) 


This  proeess  is  repeated  for  eaeh  link  in  Figure  9  using  MATLAB.  The  graphieal  result 
of  the  link  weight  ealculations  and  assignments  are  shown  in  Figure  10. 
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Figure  10.  Weighted  conflict  interference  graph  for  the  network  shown  in  Figure  8. 

Next,  we  use  the  link  weights  and  connectivity  graphs  developed  above  as  inputs 
to  Dijkstra’s  algorithm  to  calculate  the  shortest  path  based  on  interference.  Once  the 
shortest  path  has  been  determined,  we  construct  the  corresponding  routing  table  and  run 
simulations  to  compare  the  performance  of  TDMA  against  our  implementation  of 
STDMA. 


36 


VI.  SIMULATIONS,  EXPERIMENTS,  AND  PERFORMANCE 

EVALUATION 

As  discussed  in  the  preeeding  ehapter,  our  interferenee  aware  routing  method 
requires  evaluation  at  the  PHY,  DLL,  and  network  layers.  As  sueh,  we  use  two  different 
eomputing  tools  to  assist  in  the  task  of  simulating  our  system.  The  first  step  is  to 
ealeulate  SIR  values  and  develop  a  weighted  eonfliet  graph  based  on  interferenee.  This 
step  is  aeeomplished  using  MATLAB,  a  eommon  engineering  and  mathematies 
programming  language.  The  seeond  step  is  to  setup  a  WSN  and  use  least-eost  routing 
with  both  TDMA  and  STDMA  so  that  eomparisons  ean  be  made  on  the  performanee  of 
our  WSN  routing  teehnique.  For  the  seeond  set  of  tasks,  we  use  QualNet  [24],  whieh  is  a 
eomprehensive  suite  of  tools  for  modeling  wired  and  wireless  networks  ineluding  WSNs. 

A,  LEAST  COST  PATH  GENERATION  IN  MATLAB 

An  overview  of  the  steps  aeeomplished  in  MATLAB,  detailed  in  Chapter  V,  ean 
be  seen  in  Figure  1 1 .  The  speeifie  MATLAB  seripts  used  to  implement  eaeh  bloek  ean 
be  found  in  Appendix  A  through  Appendix  F. 


Figure  1 1 .  MATLAB  flowehart  of  least-eost  path  generation. 


We  simulate  two  network  eonfigurations,  one  with  10  nodes  and  another  with  21 
nodes.  The  souree  and  destination  node  were  arbitrarily  ehosen  to  be  node  1  and  node  10 
in  the  ten-node  network,  respeetively.  Similarly,  the  souree  and  destination  node  were 
ehosen  to  be  node  1  and  node  21  in  the  21 -node  network,  respeetively.  We  assume  in 
these  eases  that  nodes  1  and  21  are  the  gateways  for  their  respeetive  networks.  This 
allows  data  to  traverse  as  mueh  of  the  network  as  possible  and  allows  a  more  thorough 
and  realistie  evaluation  of  the  effeets  of  interferenee  on  network  performanee. 


37 


In  order  to  obtain  the  least-eost  routing  paths  based  on  interferenee,  we  followed 
the  MATLAB  flowehart  as  shown  in  Figure  11.  First  we  built  a  ten-node  and  21 -node 
network  and  caleulated  the  distanee  between  all  nodes  using  Eq.  (3.1).  Next,  we 
ealeulated  the  reeeived  power  at  eaeh  node  based  on  transmissions  from  every  other  node 
and  used  this  information  to  ealeulate  the  SIR  values  for  eaeh  node  per  Eqs.  (5.1)  and 
(5.2),  respeetively.  Finally,  we  used  the  SIR  information  to  ealeulate  the  link  metrie  eost 
using  Eq.  (5.3).  We  used  all  the  link  metrie  eosts  and  network  topology  as  inputs  to  the 
shortest  path  routing  funetion  of  MATLAB.  The  results  of  the  ten-node  and  21 -node 
shortest  path  determinations  are  shown  in  Figure  12  and  Figure  13. 


Network  topology 


Figure  12.  Ten-node  interference-based  shortest  paths. 
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Figure  13.  21 -node  interference  based  shortest  paths. 


Note  that  in  the  ten-node  case  the  least  cost  route  is  also  the  least  hop  route, 
unlike  the  21 -node  case  where  the  least  cost  route  has  more  hops  than  the  least  hop  route. 

B.  ROUTING  SIMULATION  IN  QUALNET 

QualNet  [29]  allows  for  scenario  based  simulation  of  various  network  types  of 
different  sizes  and  topology  configurations.  We  used  QualNet  to  establish  the  routes 
generated  by  our  MATLAB  results  and  added  in  a  variety  of  variables  to  observe  the 
effect  of  fading,  traffic  distribution  type,  and  MAC  layer  protocol  selection  on  WSN 
performance. 

QualNet  is  composed  of  the  following  tools: 

•  QualNet  Architect — A  graphical  experiment  design  and  visualization  tool. 
Architect  has  two  modes:  design  mode  for  designing  experiments,  and 
visualize  mode  for  running  and  visualizing  experiments. 

•  QualNet  Analyzer — A  graphical  statistics  analyzing  tool. 
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•  QualNet  Packet  Tracer — graphical  tool  to  display  and  analyze  packet 
traces. 

•  QualNet  File  Editor — A  text  editing  tool. 

•  QualNet  Command  Line  Interface — Command  line  access  to  the 
simulator. 

In  order  to  ensure  that  our  results  correspond  to  the  scenario  that  was  built  in 
MATLAB,  we  built  the  equivalent  scenario  in  QualNet,  which  is  depicted  in  Figure  14 
for  the  21-node  scenario.  In  Figure  14,  the  cloud  icon  and  blue  dashed  lines  represent  a 
wireless  subnet  connection  for  all  nodes  in  the  WSN.  In  addition  to  node  placement,  all 
wireless  properties  of  the  network  were  assigned  to  match  the  MATLAB  parameters  (i.e., 
transmit  power  and  path-loss  model). 
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Ligure  14.  21-node  scenario  built  in  QualNet. 

1,  Simulation  Parameters 

The  simulation  parameters  common  to  all  baseline  experiments  are  enumerated  in 
Tables  1  through  3,  for  the  PHY,  MAC,  and  Network  layers,  respectively. 
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Table  1.  PHY  layer  simulation  parameters. 


Packet  Reception  Model 

SIR  Threshold 

SIR  Threshold 

10 

Transmission  Power 

10  dBm 

Data  Rate 

2  Mbps 

Antenna  Gain 

OdB 

Antenna  Noise  Temperature 

290  K 

The  PHY  parameters  that  are  pertinent  to  the  operation  of  our  network  simulation 
are  presented  in  Table  1.  These  speeifieations  require  a  received  SIR  greater  than  10  in 
order  for  a  packet  to  be  considered  successfully  received  (note  that  the  SIR  threshold  P 
is  unitless).  We  assign  10  mW  (10  dBm)  transmission  power  for  all  nodes  and  a  data  rate 
of  i?  =  2  Mbps.  Our  antenna  parameters  are  a  0  dB  gain  and  290  K  antenna  temperature, 
which  is  the  standard  antenna  temperature  for  terrestrial  antenna  operation. 


Table  2.  MAC  layer  simulation  parameters. 


MAC  Protocol 

TDMA  or  STDMA 

Slot  Duration 

10  ms 

Guard  Time 

0  ms 

Inter-frame  Time 

1  ms 

Slots  per  frame 

9  or  20 

Our  MAC  parameters  specify  a  10  ms  slot  per  node  with  no  guard  time  between 
node  time  slots  and  a  1  ms  guard  between  successive  frames.  The  nine  or  20  time  slots 
specify  one  slot  per  node  per  frame  in  the  10-  and  21 -node  scenarios,  respectively. 


Table  3.  Network  layer  simulation  parameters. 


Routing  Protocol 

Static  Routes 

Networking  Protocol 

Internet  Protocol  version  four 

Our  network  routing  protocol  consists  of  using  the  shortest  paths  created  using  the 
Dijkstra  algorithm  in  MATLAB.  The  static  routes  specified  correspond  to  the  routes  in 
Figure  12  and  Figure  13.  All  nodes  utilize  standard  internet  protocol  addressing. 
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In  addition  to  these  general  parameters,  there  are  additional  parameters  that  must 
be  speeified  depending  on  the  scenario.  In  some  scenarios  we  include  Rayliegh  fading,  in 
which  case  the  fading  parameters  are  shown  in  Table  4. 


Table  4.  Fading  channel  parameters  used  in  the  simulations. 


Fading  Model 

Rayleigh 

Signal  Propagation  Speed 

3x  10^  m/s 

Propagation  Limit 

-111  dBm 

The  Rayleigh  fading  model  is  a  statistical  model  to  represent  the  fast  variation  of 
signal  amplitude  at  the  receiver.  In  wireless  propagation,  Rayleigh  fading  occurs  when 
there  is  no  line-of-sight  between  the  transmitter  and  receiver  [29]. 

2.  Experimental  Results 

The  results  of  the  experiments  for  various  combinations  of  the  MAC  protocol, 
fading  channel,  and  traffic  distribution  are  shown  in  Table  5  for  the  10-node  scenario. 
The  results  of  interest  are  the  delay  and  end-to-end  throughput  from  source  to  destination 
for  constant  bit  rate  (CBR)  and  variable  bit  rate  (VBR)  traffic  applications.  The  VBR 
traffic  model  uses  an  exponential  distribution.  The  last  two  columns  of  Table  5  indicate 
the  percent  improvement  from  using  STDMA  versus  TDMA  for  the  throughput  and  delay 
results  shown. 
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Table  5.  10-node  network.  Pereent  improvement  in  delay  and  throughput  for  the  10-node 
network  when  using  STDMA  versus  STDMA  for  varying  seenarios. 


MAC 

Protocol 

Fading  Model 

Traffic 

Model 

Throughput 

(kbps) 

Delay 

(sec) 

Throughput 

Improvement 

using 

STDMA  (%) 

Delay 

Improvement 

using 

STDMA 

(%) 

TDMA 

None 

CBR 

4283 

0.1237 

-0.4 

24.3 

STDMA 

None 

CBR 

4266 

0.0937 

TDMA 

Rayleigh 

CBR 

4105 

0.1194 

3.9 

21.5 

STDMA 

Rayleigh 

CBR 

4266 

0.0937 

TDMA 

None 

VBR 

5936 

0.1493 

3.4 

31.0 

STDMA 

None 

VBR 

5733 

0.1030 

TDMA 

Rayleigh 

VBR 

5576 

0.1488 

-0.004 

29.8 

STDMA 

Rayleigh 

VBR 

5554 

0.1044 

The  same  simulation  eombinations  were  repeated  for  the  21 -node  network  and  the 
results  are  shown  in  Table  6. 


Table  6.  Results  eomparing  delay  and  throughput  when  using  STDMA  versus  TDMA 
under  different  fading  and  traffie  models  for  a  21 -node  network. 


MAC 

Protocol 

Fading  Model 

Traffic 

Model 

Throughput 

(kbps) 

Delay 

(sec) 

Throughput 

Improvement 

using 

STDMA  (%) 

Delay 

Improvement 

using 

STDMA 

(%) 

TDMA 

None 

CBR 

4388 

0.2057 

-6.7 

25.9 

STDMA 

None 

CBR 

4095 

0.1524 

TDMA 

Rayleigh 

CBR 

4388 

0.2057 

-10.7 

25.9 

STDMA 

Rayleigh 

CBR 

3917 

0.1524 

TDMA 

None 

VBR 

6088 

0.3919 

-20.3 

53.1 

STDMA 

None 

VBR 

4850 

0.1839 

TDMA 

Rayleigh 

VBR 

6088 

0.3919 

-20.3 

53.1 

STDMA 

Rayleigh 

VBR 

4850 

0.1836 

In  general,  the  performanee  trends  showed  no  improvement  in  throughput  and,  in 


some  eases,  a  signilieant  deerease  in  throughput  for  the  STDMA  scenario  as  compared 
with  the  TDMA  scenario.  This  is  particularly  true  for  the  21-node  case  (Table  6). 
However,  in  all  cases  there  were  significant  decreases  in  end-to-end  delay,  which 
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indicates  there  is  value  in  the  approach  taken.  Our  interpretation  of  the  results  is  that  slot 
re-use  allows  slots  to  matriculate  faster  across  the  network  but  at  the  expense  of 
interference  causing  the  received  SIR  to  drop  below  the  required  threshold.  This  loss  of 
packets  results  in  a  decreased  throughput. 

In  an  attempt  to  improve  the  throughput  performance,  we  re-attempted  the  21- 
node  CBR  no  fading  scenario  with  varying  values  of  SIR  threshold.  We  simulated  two 
sets  of  oases:  (1)  small  changes  around  our  initial  threshold  value  of  10  and  (2)  threshold 
ohanges  over  a  larger  scale.  The  results  when  the  SIR  is  varied  on  a  small  scale  are 
shown  in  Table  7,  whereas  the  results  when  the  SIR  is  varied  on  a  large  scale  are  shown 
in  Table  8.  The  values  in  the  Tables  7  and  8  indioate  a  decreasing  trend  in  throughput  as 
the  SIR  threshold  becomes  more  discriminatory.  In  addition,  we  did  not  observe  any 
inorease  in  throughput  as  compared  to  our  throughput  for  the  STDMA  case  with  no 
fading  and  CBR  traffic  when  the  SIR  threshold  was  lowered  below  ten.  These  results 
confirm  that  the  STDMA  approaeh  is  not  beneficial  in  terms  of  throughput  but  is  useful 
in  reducing  transmission  delay. 


Table  7.  Small  variations  in  SIR  threshold  for  throughput  analysis  of  the  21 -node  network. 


SIR  Threshold 

Throughput  (kbps) 

1 

4095 

2 

4095 

3 

4095 

4 

4095 

5 

4095 

6 

3917 

7 

4095 

8 

3917 

9 

3739 

10 

3917 

11 

3739 

12 

3739 

13 

3917 

14 

3739 

15 

3561 
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Table  8. 


Large  variations  in  SIR  Threshold  for  throughput  analysis  of  the  21 -node 

network. 


SIR  Threshold 

Throughput  (kbps) 

5 

4095 

10 

3917 

15 

3561 

20 

3561 

25 

3561 

30 

3383 

35 

2493 

40 

1638 

45 

0 

50 

0 

In  this  chapter  we  implemented  the  interference  aware  STDMA  and  least-cost 
routing  method  proposed  in  Chapter  V.  The  interference  introduced  by  the  transmission 
groups  caused  poor  network  throughput  performance,  particularly  in  the  21 -node 
network.  Additionally,  we  performed  simulations  to  determine  if  we  could  gain  a 
throughput  increase  by  manipulating  the  SIR  threshold  but  did  not  see  any  improvement 
for  less  discriminant  choices  of  SIR  threshold.  We  saw  a  sizeable  improvement  in 
network  delay  by  utilizing  transmission  group  scheduling  and  interference  aware  least 
cost  routing. 
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VII.  CONCLUSIONS 


A.  CONCLUDING  REMARKS 

Throughout  the  literature  there  is  ample  proof  of  the  value  of  using  a  STDMA 
MAC  protoeol  for  improved  network  performanee.  In  a  WSN,  this  ean  be  exploited  to 
better  utilize  preeious  and  limited  network  resourees  and  inerease  the  operating  network 
lifetime;  however,  many  STDMA  implementations  incur  a  costly  computational  penalty, 
which  can  mitigate  the  savings  benefits  of  the  optimized  schedules  they  provide. 

In  this  thesis,  we  presented  the  idea  of  a  simplified  scheduling  philosophy  using 
STDMA.  As  interference  is  a  bottleneck  to  network  performance,  we  studied  the  effect  of 
interference  when  using  STDMA  integrated  with  a  PhIM.  The  implementation  of  the 
PhIM  allowed  us  to  capture  link  interference  metrics.  We  used  the  concept  of  spatial 
reuse  to  schedule  simultaneous  transmissions  by  implementing  a  conflict  graph.  We 
implemented  Dijkstra’s  algorithm,  where  the  link  costs  are  obtained  from  SIR 
calculations  determined  from  the  conflict  graph.  Using  these  links,  we  determined  the 
shortest  path  from  a  sensor  node  to  the  gateway.  Via  our  simulations  using  MATLAB 
and  QualNet,  we  showed  that  our  proof-of-concept  improves  delay  in  the  network.  An 
advantage  to  our  proposed  scheme  is  that  it  does  not  require  heuristic  algorithms  or 
messaging  overhead  for  proper  system  operation;  however,  a  tradeoff  is  that  our  method 
did  not  help  improve  throughput.  Nonetheless,  we  see  the  results  obtained  as  a  proof-of- 
concept  that  a  deterministic  PhIM  can  provide  value  in  network  design  and  performance 
despite  its  algorithmic  computational  complexity. 

B,  FUTURE  WORK 

1,  Network  Generation  Automation  and  QualNet  Integration 

In  order  to  create  a  more  flexible  and  robust  scenario  generation  environment,  we 
suggest  the  implementation  of  automated  network  generation.  This  should  include 
randomized  node  placement  and  calculation  of  power  and  interference  values  for 
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arbitrary  node  distances.  This  generalized  approach  would  permit  more  interesting  and 
less  restrictive  scenario  development  and  testing. 

The  QualNet  simulation  environment  provides  a  comprehensive  simulation 
environment  that  is  designed  to  be  customizable  and  allow  high-level  development  at  all 
network  layers.  In  order  to  establish  a  more  generalized  network  topology  and  integrate 
the  PhIM  functionality,  significant  modification  to  the  PHY  layer  functions  included  in 
QualNet  is  required.  By  accomplishing  the  necessary  integration  steps,  very  large 
scenarios  can  be  developed  with  varying  node  density  to  obtain  data  to  verify  the 
usefulness  of  our  approach. 

2,  Development  of  a  Formal  Protocol 

In  order  for  a  network  to  operate,  there  must  be  a  set  of  rules  for  the  constituent 
devices  to  operate  under.  This  requires  the  formalization  of  inter-layer  communications 
to  pass  SIR  information  from  the  PHY  layer  to  the  network  layer  to  be  used  for  routing. 
Additionally,  the  development  of  messaging  components  including  message  header  and 
internal  packet  structures  are  needed  to  disseminate  information  for  network  scheduling 
among  the  nodes. 

3.  Testing  in  a  Real-World  WSN 

Finally,  building  and  testing  a  WSN  using  the  scheduling  considerations  of  the 
thesis  and  verifying  delay  improvement  shown  in  the  simulation  results  obtained  in 
QualNet  would  provide  practical  analysis  to  our  proof-of-concept  developed  in  this 
thesis.  The  information  gained  by  a  real-world  test  can  be  used  to  determine  if  the 
interference  calculations  are  overly  optimistic  or  pessimistic  and  perhaps  lead  to  a 
relaxation  of  some  of  our  simplifying  assumptions. 
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APPENDIX 


The  appendices  that  follow  document  the  MATLAB  code  utilized  for  creation  of 
the  interference  based  least  cost  network  path. 

A.  NODE  TOPOLOGY  FUNCTION  FOR  lO-NODE  NETWORK 

This  function  serves  to  populate  a  ten-node  network  with  all  nodes  separated  by 

10  meters.  The  maxx  and  maxy  parameters  are  used  to  scale  the  distance  and  the  node 

array  is  used  for  manual  entry  of  node  coordinates.  There  is  a  drawFigure  argument  that 

can  be  used  to  produce  a  graphical  depiction  of  the  network  topology. 

function  [node]  =  topolO (maxx,  maxy,  drawFigure) ; 

%  Generate  network  topology 

%  Maxx  =  maximum  x  distance  for  a  node  to  be  placed 
%  Maxy  =  maximum  y  distance  for  a  node  to  be  placed 
%  drawFigure  =  boolean  input  1  for  plot  0  for  no  plot 

node  =  [1  1;2  1;3  1;1  2;2  2 ; 3  2 ; 1  3;2  3;3  3;2  4], *0.1;  %Enter 
Coordinates  of  nodes  here 
node ( : , 1 )  =  node ( : , 1 ) *maxx; 
node(:,2)  =  node ( : , 2 ) *maxy; 


if  drawFigure  >=  1 

colordef  none,  whitebg 

figure  ( 1 ) ; 

axis  equal 

hold  on; 

box  on; 

plot (node (: ,  1),  node ( : ,  2),  'k.',  'MarkerSize' ,  5); 

title (  'Network  topology' ) ; 

xlabel ( 'X  meters' ) ; 

ylabel ( 'Y  meters' ) ; 

axis  ([0,  40,  0,  50]); 

set(gca,  'XTick' ,  [0;  40]); 

set(gca,  'YTick' ,  [50]); 

end 

return; 


B,  NODE  TOPOLOGY  FUNCTION  FOR  2 1-NODE  NETWORK 

This  function  serves  to  populate  a  21 -node  network  with  all  nodes  separated  by  10 
meters.  The  maxx  and  maxy  parameters  are  used  to  scale  the  distance  and  the  node  array 
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is  used  for  manual  entry  of  node  coordinates.  There  is  a  drawFigure  argument  that  can  be 
used  to  produce  a  graphical  depiction  of  the  network  topology. 

function  [node]  =  topo2 1 (maxx,  maxy,  drawFigure); 

%  Generate  network  topology 

%  Maxx  =  maximum  x  distance  for  a  node  to  be  placed 
%  Maxy  =  maximum  y  distance  for  a  node  to  be  placed 
%  drawFigure  =  boolean  input  1  for  plot  0  for  no  plot 

node  =  [1  1;2  1;3  1;4  1;5  1;1  2;2  2 ; 3  2 ; 4  2 ; 5  2 ; 1  3;2  3;3  3;4  3;5  3;1 

4;  .  .  . 

2  4 ; 3  4 ; 4  4 ; 5  4 ; 3  5]. *0.1;  %Enter  Coordinates  of  nodes  here 
node ( : , 1 )  =  node ( : , 1 ) *maxx; 
node(:,2)  =  node ( : , 2 ) *maxy; 

if  drawFigure  >=  1 

colordef  none,  whitebg 

figure  ( 1 ) ; 

axis  equal 

hold  on; 

box  on; 

plot (node (: ,  1),  node ( : ,  2),  'k.',  'MarkerSize' ,  5); 

title (  'Network  topology' ) ; 

xlabel ( 'X  meters' ) ; 

ylabel (  'Y  meters' ) ; 

axis  ([0,  60,  0,  60]); 

setigca,  'XTick' ,  [0;  60]); 

set(gca,  'YTick' ,  [60]); 

end 

return; 


C.  DISTANCE  CALCULATION  FUNCTION 

This  function  is  called  to  calculate  the  Euclidean  distance  between  nodes. 

function  d  =  node_dist (xl , yl , x2  ,  y2 )  ; 

%determines  the  distance  beween  nodes 
d  =  sqrt  (  (x2-xl ) ''2+ (y2-yl ) ''2  )  ; 
return; 

D,  RECEIVED  POWER  CALCULATION  FUNCTION 

This  function  is  called  to  calculate  the  received  power  based  on  distance  between 
nodes  using  the  path-loss  model. 

function  Prx  =  pathlosscalc ( Ptx, d) ; 

%Determines  change  in  power  due  to  distance  losses 

Prx  =  Ptx*d''(-2);  %2  is  the  path  loss  factor  which  is  typical 

%Prx  =  loglO(Prx) 

return 
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E. 


MATLAB  SCRIPT  TO  DETERMINE  THE  LEAST  COST  PATH  FROM 
NODE  1  TO  NODE  10  IN  THE  TEN-NODE  NETWORK. 


This  script  calculates  the  SIR  and  link  values  as  well  as  the  eonneetivity  and 
conflict  graphs  for  the  ten-node  network.  The  function  then  uses  the  built  in 
shortestpathO  function  in  MATLAB  to  determine  the  shortest  path  and  associated  path 
cost  through  the  network. 

%5  node  SINK  link-weight  generation 

clear 

%setup  constansts 

Ptx  =  10*10^' (-3)  ;  %Transmit  Power  of  lOmW 


%Intiate  Topology 

nodes  =  topolO ( 100 , 100 , 0 ) ;  %returns  coordinates  of  nodes  in  order 
%calculate  distances  between  nodes 


xl  = 

nodes 

(1, 

1) 

yl  = 

nodes 

(1, 

2) 

x2  = 

nodes 

(2, 

1) 

y2  = 

nodes 

(2, 

2) 

x3  = 

nodes 

(3, 

1) 

y3  = 

nodes 

(3, 

2) 

x4  = 

nodes 

(4, 

1) 

y4  = 

nodes 

(4, 

2) 

x5  = 

nodes 

(5, 

1) 

y5  = 

nodes 

(5, 

2) 

x6  = 

nodes 

(6, 

1) 

y6  = 

nodes 

(6, 

2) 

x7  = 

nodes 

(7, 

1) 

y7  = 

nodes 

(7, 

2) 

x8  = 

nodes 

(8, 

1) 

y8  = 

nodes 

(8, 

2) 

x9  = 

nodes 

(9, 

1) 

y9  = 

nodes 

(9, 

2) 

xlO  = 

=  nodes 

(1 

0, 

ylO  = 

=  nodes 

(1 

0, 

%distance  from  each  node  to  each  other  node 

dl2  =  node_dist (xl , yl , x2  ,  y2 )  ; 
dl3  =  node_dist (xl , yl , x3 ,  y3 )  ; 
dl4  =  node_dist (xl , yl , x4 , y4 ) ; 
dl5  =  node_dist (xl , yl ,  x5 ,  y5 )  ; 
dl6  =  node_dist (xl , yl ,  x6 ,  y6 )  ; 
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dl7  =  node_dist (xl , yl , x7 , y7 ) ; 
dl8  =  node_dist (xl , yl ,  x8 ,  y8 )  ; 
dl9  =  node_dist (xl ,  yl ,  x9 ,  y9 )  ; 
dllO  =  node_dist (xl , yl ,  xlO ,  ylO )  ; 

d21  =  dl2; 

d23  =  node_dist (x2 , y2 , x3 , y3 ) ; 
d24  =  node_dist (x2 , y2 , x4 , y4 ) ; 
d25  =  node_dist (x2 ,  y2 ,  x5 ,  y5 )  ; 
d2  6  =  node_dist (x2 ,  y2 ,  x6 ,  y6 )  ; 
d27  =  node_dist (x2 , y2 , x7 , y7 ) ; 
d2  8  =  node_dist (x2 ,  y2 ,  x8 ,  y8 )  ; 
d2  9  =  node_dist (x2 ,  y2 ,  x9 ,  y9 )  ; 
d210  =  node_dist (x2 , y2 ,  xlO ,  ylO )  ; 

d31  =  dl3; 
d32  =  d23; 

d34  =  node_dist (x3 , y3 , x4 , y4 ) ; 
d35  =  node_dist (x3 , y3 ,  x5 ,  y5 )  ; 
d36  =  node_dist (x3 ,  y3 ,  x6 ,  y6 )  ; 
d37  =  node_dist (x3 , y3 , x7 ,  y7 ) ; 
d38  =  node_dist (x3 , y3 ,  x8 ,  y8 )  ; 
d3  9  =  node_dist (x3 ,  y9 ,  x5 ,  y9 )  ; 
d310  =  node_dist (x3, y3, xlO, ylO) ; 

d41  =  dl4; 
d42  =  d24; 
d43  =  d34; 

d45  =  node_dist (x4 , y4 , x5 , y5 ) ; 
d46  =  node_dist (x4 , y4 , x6 , y6 ) ; 
d47  =  node_dist (x4  ,  y4  ,  x7  ,  y7  )  ; 
d4  8  =  node_dist (x4 ,  y4 ,  x8 ,  y8 )  ; 
d4  9  =  node_dist (x4 , y4 ,  x9 ,  y9 )  ; 
d410  =  node_dist (x4 , y4 ,  xlO ,  ylO )  ; 

d51  =  dl5; 
d52  =  d25; 
d53  =  d35; 
d54  =  d45; 

d56  =  node_dist (x5 , y5 , x6 , y6 ) ; 
d57  =  node_dist (x5 , y5 ,  x7 ,  y7 )  ; 
d58  =  node_dist (x5 ,  y5 ,  x8 ,  y8 )  ; 
d59  =  node_dist (x5 , y5 ,  x9 ,  y9 )  ; 
dSlO  =  node_dist (x5, y5, xlO, ylO) ; 

d61  =  dl6; 
d62  =  d26; 
d63  =  d36; 
d64  =  d46; 
d65  =  d56; 

d67  =  node_dist (x6 , y6 , x7 , y7 ) ; 
d68  =  node_dist (x6 , y6 , x8 ,  y8 )  ; 
d69  =  node_dist (x6 ,  y6 ,  x9 ,  y9 )  ; 
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d610  = 

=  node  dist (x6 , y6 , xl 

0,  ylO 

d71  = 

dl7; 

d72  = 

d2  7; 

d73  = 

d37; 

d74  = 

d4  7; 

d75  = 

d57; 

d7  6  = 

d67; 

d7  8  = 

node  dist (x7 , y7 , x8 , 

y8)  ; 

d7  9  = 

node  dist (x7 , y7 , x9 , 

y9)  ; 

d710  = 

=  node  dist (x7 , y7 , xl 

0,  ylO 

d81  = 

dl8; 

d82  = 

d2  8; 

d83  = 

d38; 

d84  = 

d4  8; 

d85  = 

d58; 

d8  6  = 

d68; 

d87  = 

d7  8; 

d8  9  = 

node  dist (x8 , y8 ,  x9 , 

y9)  ; 

d810  = 

=  node  dist (x8 , y8 ,  xl 

0,  ylO 

d91  = 

dl9; 

d92  = 

d2  9; 

d93  = 

d39; 

d94  = 

d4  9; 

d95  = 

d59; 

d96  = 

d69; 

d97  = 

d7  9; 

d98  = 

d8  9; 

d910  = 

=  node  dist (x9 , y9 , xl 

0,  ylO 

dlOl  = 

=  dllO; 

dl02  = 

=  d210; 

dl03  = 

=  d310; 

dl04  = 

=  d410; 

dl05  = 

=  dSlO; 

dl06  = 

=  d610; 

dl07  = 

=  d710; 

dl08  = 

=  d810; 

dl09  = 

=  d910; 

%Calculate  the  recieved  power 

Prxl2 

=  pathlosscalc ( Ptx, 

dl2)  ; 

Prxl3 

=  pathlosscalc ( Ptx, 

dl3)  ; 

Prxl4 

=  pathlosscalc ( Ptx, 

dl4)  ; 

PrxlS 

=  pathlosscalc ( Ptx, 

dl5)  ; 

Prxl  6 

=  pathlosscalc ( Ptx, 

dl6)  ; 

Prxl7 

=  pathlosscalc ( Ptx, 

dl7)  ; 

Prxl  8 

=  pathlosscalc ( Ptx, 

dl8)  ; 

Prxl  9 

=  pathlosscalc ( Ptx, 

dl9)  ; 

PrxllC 

1  =  pathlosscalc ( Ptx 

;,dll0 
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Prx21  =  pathlosscalc ( Ptx, d2 1 ) ; 
Prx23  =  pathlosscalc (Ptx, d23)  ; 
Prx24  =  pathlosscalc (Ptx, d24) ; 
Prx25  =  pathlosscalc (Ptx, d25) ; 
Prx26  =  pathlosscalc ( Ptx, d2 6 ) ; 
Prx27  =  pathlosscalc (Ptx, d27) ; 
Prx28  =  pathlosscalc ( Ptx, d2 8 ) ; 
Prx29  =  pathlosscalc ( Ptx, d2 9 ) ; 
Prx210  =  pathlosscalc (Ptx, d210) ; 

Prx31  =  pathlosscalc (Ptx, d31) ; 
Prx32  =  pathlosscalc ( Ptx, d32 ) ; 
Prx34  =  pathlosscalc ( Ptx, d34 ) ; 
Prx35  =  pathlosscalc ( Ptx, d35 ) ; 
Prx36  =  pathlosscalc (Ptx, d36) ; 
Prx37  =  pathlosscalc ( Ptx, d37 ) ; 
Prx38  =  pathlosscalc (Ptx, d38) ; 
Prx39  =  pathlosscalc (Ptx, d39) ; 
Prx310  =  pathlosscalc (Ptx, d310) ; 

Prx41  =  pathlosscalc ( Ptx, d4 1 ) ; 
Prx42  =  pathlosscalc ( Ptx, d42 ) ; 
Prx43  =  pathlosscalc (Ptx, d43)  ; 
Prx45  =  pathlosscalc (Ptx, d45) ; 
Prx46  =  pathlosscalc ( Ptx, d4 6 ) ; 
Prx47  =  pathlosscalc (Ptx, d47) ; 
Prx48  =  pathlosscalc ( Ptx, d4 8 ) ; 
Prx49  =  pathlosscalc ( Ptx, d4 9 ) ; 
Prx410  =  pathlosscalc (Ptx,  d410)  ; 


PrxSl  =  pathlosscalc (Ptx, dSl) ; 
Prx52  =  pathlosscalc ( Ptx, d52 ) ; 
Prx53  =  pathlosscalc ( Ptx, d53 ) ; 
Prx54  =  pathlosscalc ( Ptx, d54 ) ; 
Prx56  =  pathlosscalc (Ptx, d56) ; 
Prx57  =  pathlosscalc ( Ptx, d57 )  ; 
Prx58  =  pathlosscalc (Ptx, d58) ; 
Prx59  =  pathlosscalc (Ptx, d59) ; 
PrxSlO  =  pathlosscalc (Ptx, dSlO) ; 

PrxGl  =  PrxlG; 

Prx62  =  Prx26; 

Prx63  =  Prx36; 

Prx64  =  Prx46; 

Prx65  =  Prx56; 

Prx67  =  pathlosscalc ( Ptx, d67 ) ; 
Prx68  =  pathlosscalc (Ptx,  d68)  ; 
Prx69  =  pathlosscalc (Ptx, d69) ; 
PrxGlO  =  pathlosscalc (Ptx, d610) ; 

Prx71  =  Prxl7; 

Prx72  =  Prx27; 


Prx73  = 

Prx37 ; 

Prx74  = 

Prx47; 

Prx75  = 

Prx57 ; 

Prx76  = 

Prx67 ; 

Prx78  = 

pathlosscalc ( Ptx, d7  8 ) ; 

Prx79  = 

pathlosscalc( Ptx, d7  9 ) ; 

Prx710 

=  pathlosscalc (Ptx, d710) ; 

Prx81  = 

Prxl8; 

Prx82  = 

Prx2  8 ; 

Prx83  = 

Prx38 ; 

Prx84  = 

Prx48; 

Prx85  = 

Prx58 ; 

Prx86  = 

Prx68 ; 

Prx87  = 

Prx7  8 ; 

Prx89  = 

pathlosscalc ( Ptx,  d8 9 )  ; 

Prx810 

=  pathlosscalc (Ptx, d810) ; 

Prx91  = 

Prxl9; 

Prx92  = 

Prx2  9 ; 

Prx93  = 

Prx39; 

Prx94  = 

Prx4  9 ; 

Prx95  = 

Prx59; 

Prx96  = 

Prx69; 

Prx97  = 

Prx7  9 ; 

Prx98  = 

Prx89; 

Prx910 

=  pathlosscalc (Ptx, d910) ; 

PrxlOl 

=  Prxl 10; 

Prxl02 

=  Prx210; 

Prxl03 

=  Prx310; 

Prxl04 

=  Prx410; 

PrxlOS 

=  PrxSlO; 

PrxlOG 

=  PrxGlO; 

Prxl07 

=  Prx710; 

Prxl08 

=  Prx810; 

Prxl09 

=  Prx910; 

%Calculate  Signal  to  Interference  Ratios 

%SIR  for  transmissions  from  1  to  2,4,  or  5 

SIR2froml  =  ( Prx2 1 ) / ( Prx32  +  Prx72  +  Prx92 ) ; 
SIR4froml  =  (Prx41) / (Prx34+Prx74+Prx94) ; 
SIRSfroml  =  (PrxSl) / (Prx35+Prx75+Prx95) ; 


%SIR  for  transmissions  from  2  to  1,3, 4, 5, 6 

SIRlfrom2  =  Prxl2 /Prxl 8 ; 

SIR3from2  =  Prx32/Prx38; 

SIR4from2  =  Prx42 /Prx4 8 ; 
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SIR5from2  =  Prx52/Prx58; 
SIR6from2  =  Prx62/Prx68; 


%SIR  for  transmissions  from  3  to  2,5,6 

SIR2from3  =  Prx23/ (Prx21+Prx27+Prx29) ; 
SIR5from3  =  Prx53/ (Prx51+Prx57+Prx59) ; 
SIR6from3  =  Prx63/ (Prx61+Prx67+Prx69) ; 


%SIR  for  transmissions  from  4 


to  1 , 2 , 5 , 7 , 8 


SIRlfrom4 
SIR2f rom4 
SIR5f rom4 
SIR7f rom4 
SIR8from4 


Prx41/Prx61; 
Prx42 / Prx62 ; 
Prx45/Prx65; 
Prx47/Prx67; 
Prx48/Prx68; 


%SIR  for  transmissions 


%Note,  5  ha 
SIRlfromS  = 
SIR2from5  = 
SIR3from5  = 
SIR4from5  = 
SIR6from5  = 
SIR7from5  = 
SIR8from5  = 
SIR9from5  = 


no  two  hop 
Prx51 ; 

Prx52 ; 
Prx53; 

Prx54 ; 
Prx56; 

Prx57 ; 

Prx58 ; 
Prx59; 


from  5  to  1 , 2 , 3 , 4 , 6 , 7 , 8 , 9 
neighbors 


%SIR  for  transmissions  from  6  to  2, 3, 5, 8, 9 
SIR2from6  =  Prx62 /Prx42 ; 

SIR3from6  =  Prx63/Prx43; 

SIR5from6  =  Prx65/Prx45; 

SIR8from6  =  Prx68/Prx48; 

SIR9from6  =  Prx69/Prx49; 

%SIR  for  transmissions  from  7  to  4,5,8,10 
SIR4from7  =  Prx74/ (Prxl4+Prx34+Prx94) ; 
SIR5from7  =  Prx75/ (Prxl5tPrx35tPrx95) ; 
SIR8from7  =  Prx87/ (Prxl8+Prx38+Prx98) ; 
SIR10from7  =  Prx710/ (Prxll0tPrx310tPrx910) ; 

%SIR  for  transmissions  from  8  to  4,5,6,7,9,10 
SIR4from8  =  Prx84/Prx24; 

SIR5from8  =  Prx85/Prx25; 

SIR6from8  =  Prx86/Prx26; 

SIR7from8  =  Prx87/Prx27; 

SIR9from8  =  Prx89/Prx29; 

SIR10from8  =  Prx810/Prx210; 

%SIR  for  transmissions  from  9  to  5,6,8,10 
SIR5from9  =  Prx95/ (Prxl5tPrx35+Prx75) ; 
SIR6from9  =  Prx96/ (Prxl6+Prx36+Prx76) ; 
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SIR8from9  =  Prx98/ (Prxl8+Prx38+Prx78) ; 

SIR10from9  =  Prx910/ (Prxll0+Prx310+Prx710) ; 

%Fill  out  the  Link  Weight  matrix  for  shortest  path  input 

%Note  nodes  beyond  one  hop  are  considered  disconnected  and  assigned  a 

link 

%weight  of 


W  =  ([SIR4froml  SIR2froml  SIRSfroml  SIR4from2  SIR5from2  SIR6from2  ... 
SIR2from3  SIR5from3  SIR6from3  SIR5from4  SIR5from6  SIR8from4  . . . 
SIR8from6  SIR7from4  SIR9from6  SIR8from7  SIR8from9  SIR10from7  ... 
SIR10from8  SIR10from9 

DG  =  sparse ([1  1122233346464679789  10],  [4  254562 
5  6  5  5  8  8  7  9  8  8  10  10  10  9] ,W) ; 

%h  =  view (biograph (DG, [],' ShowWeights' on' ) ) 

[dist,  path,  pred]  =  graphshortestpath (DG,  1,  10) 

F.  MATLAB  SCRIPT  TO  DETERMINE  THE  LEAST  COST  PATH  FROM 
NODE  1  TO  NODE  21  IN  THE  21-NODE  NETWORK. 

This  script  calculates  the  SIR  and  link  values  as  well  as  the  eonneetivity  and 
conflict  graphs  for  the  21 -node  network.  The  function  then  uses  the  built  in 
shortestpathO  function  in  MATLAB  to  determine  the  shortest  path  and  associated  path 
cost  through  the  network. 

%5  node  SINR  link-weight  generation 

clear 

%setup  constansts 

Ptx  =  10*10^ (-3);  %Transmit  Power  of  lOmW 


%Intiate  Topology 

nodes  =  topo21 (100,  100, 0)  ;  %returns  coordinates  of  nodes  in  order 


%calculate  distances  between  nodes 


xl  =  nodes  (1,1); 
yl  =  nodes  (1,2)  ; 
x2  =  nodes  (2,1); 
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y2  = 

nodes 

(2,2)  ; 

x3  = 

nodes 

(3,1)  ; 

y3  = 

nodes 

(3,2)  ; 

x4  = 

nodes 

(4,1)  ; 

y4  = 

nodes 

(4,2)  ; 

x5  = 

nodes 

(5,1)  ; 

y5  = 

nodes 

(5,2)  ; 

x6  = 

nodes 

(6,1)  ; 

y6  = 

nodes 

(6,2)  ; 

x7  = 

nodes 

(7,1)  ; 

y7  = 

nodes 

(7,2)  ; 

x8  = 

nodes 

(8,1); 

y8  = 

nodes 

(8,2); 

x9  = 

nodes 

(9,1)  ; 

y9  = 

nodes 

(9,2)  ; 

xlO 

=  nodes 

(10,1) 

ylO 

=  nodes 

(10,2) 

xll 

=  nodes 

(11,1) 

yll 

=  nodes 

(11,2) 

xl2 

=  nodes 

(12,1) 

yl2 

=  nodes 

(12,2) 

xl3 

=  nodes 

(13,1) 

yl3 

=  nodes 

(13,2) 

xl4 

=  nodes 

(14,1) 

yl4 

=  nodes 

(14,2) 

xl5 

=  nodes 

(15,1) 

yl5 

=  nodes 

(15,2) 

xl  6 

=  nodes 

(16,1) 

yl  6 

=  nodes 

(16,2) 

xl7 

=  nodes 

(17,1) 

yl7 

=  nodes 

(17,2) 

xl8 

=  nodes 

(18,1) 

yl8 

=  nodes 

(18,2) 

xl9 

=  nodes 

(19,1) 

yl9 

=  nodes 

(19,2) 

x2  0 

=  nodes 

(20,1) 

y2  0 

=  nodes 

(20,2) 

x21 

=  nodes 

(21,1) 

y21 

=  nodes 

(21,2) 

%distance  from  each  node  to  each  other  node 

dl2  =  node_dist (xl , yl , x2  ,  y2 ) ; 
dl3  =  node_dist (xl , yl , x3 ,  y3 )  ; 
dl4  =  node_dist (xl , yl , x4  ,  y4 )  ; 
dl5  =  node_dist (xl , yl ,  x5 ,  y5 )  ; 
dl6  =  node_dist (xl , yl , x6 ,  y6 ) ; 
dl7  =  node_dist (xl , yl , x7 , y7 ) ; 
dl8  =  node_dist (xl , yl ,  x8 ,  y8 )  ; 
dl9  =  node_dist (xl , yl ,  x9 ,  y9 )  ; 
dllO  =  node_dist (xl , yl , xlO ,  ylO )  ; 
dill  =  node_dist (xl , yl ,  xl 1 ,  yl 1 )  ; 
dll2  =  node_dist (xl , yl , xl2  ,  yl2 )  ; 
dll3  =  node_dist (xl, yl, xl3, yl3) ; 


58 


dll4 
dll5 
dll  6 
dll7 
dll  8 
dll  9 
dl20 
dl21 

d21  = 

d23  = 

d24  = 

d25  = 

d2  6  = 

d27  = 

d2  8  = 

d2  9  = 

d210 

d211 

d212 

d213 

d214 

d215 

d216 

d217 

d218 

d219 

d220 

d221 


d31  = 

d32  = 

d34  = 

d35  = 

d36  = 

d37  = 

d38  = 

d39  = 

d310 

d311 

d312 

d313 

d314 

d315 

d316 

d317 

d318 

d319 

d320 

d321 

d41  = 


=  node_dist (xl , yl ,  xl4  ,  yl4  )  ; 

=  node_dist (xl, yl, xl5, yl5) ; 

=  node_dist (xl ,  yl ,  xl 6 ,  yl 6 )  ; 

=  node_dist (xl , yl ,  xl7  ,  yl7  )  ; 

=  node_dist (xl , yl , xl 8 ,  yl 8 )  ; 

=  node_dist (xl , yl ,  xl 9 ,  yl 9 )  ; 

=  node_dist (xl , yl ,  x20 ,  y20 )  ; 

=  node_dist (xl , yl ,  x2 1 ,  y2 1 )  ; 

dl2; 

node_dist (x2 , y2 ,  x3 ,  y3 )  ; 
node_dist (x2 , y2 , x4  ,  y4 ) ; 
node_dist (x2 , y2 ,  x5 ,  y5 )  ; 
node_dist (x2 , y2 ,  x6 ,  y6 )  ; 
node_dist (x2 , y2 , x7  ,  y7 ) ; 
node_dist (x2 , y2 ,  x8 ,  y8 )  ; 
node_dist (x2 , y2 ,  x9 ,  y9 )  ; 

=  node_dist (x2 , y2 , xlO ,  ylO )  ; 
=  node_dist (x2 , y2  ,  xl 1 ,  yl 1 )  ; 
=  node_dist (x2 , y2 , xl2  ,  yl2 )  ; 
=  node_dist (x2, y2, xl3, yl3) ; 
=  node_dist (x2  ,  y2  ,  xl4  ,  yl4  )  ; 
=  node_dist (x2, y2, xl5, yl5) ; 
=  node_dist (x2 ,  y2 ,  xl 6 ,  yl 6 )  ; 
=  node_dist (x2 , y2 , xl7  ,  yl7 )  ; 
=  node_dist (x2 , y2 ,  xl 8 ,  yl 8 )  ; 
=  node_dist (x2 , y2 ,  xl 9 ,  yl 9 )  ; 
=  node_dist (x2 , y2 , x20 ,  y20 )  ; 
=  node_dist (x2 , y2 , x2 1 ,  y2 1 )  ; 


dl3; 
d2  3; 

node_dist (x3 , y3 ,  x4  ,  y4  )  ; 
node_dist (x3 , y3 ,  x5 ,  y5 )  ; 
node_dist (x3 , y3 ,  x6 ,  y6 )  ; 
node_dist (x3 , y3 ,  x7 ,  y7 )  ; 
node_dist (x3 , y3 ,  x8 ,  y8 )  ; 
node_dist (x3,y9,x5,y9)  ; 

=  node_dist (x3, y3, xlO, ylO) ; 
=  node_dist (x3 , y3 ,  xl 1 ,  yl 1 )  ; 
=  node_dist (x3 , y3 , xl2  ,  yl2  )  ; 
=  node_dist (x3, y3, xl3, yl3) ; 
=  node_dist (x3, y3, xl4, yl4) ; 
=  node_dist (x3, y3,  xl5,  yl5)  ; 
=  node_dist (x3 , y3 ,  xl 6 ,  yl 6 )  ; 
=  node_dist (x3, y3, xl7, yl7) ; 
=  node_dist (x3 , y3 ,  xl 8 ,  yl 8 )  ; 
=  node_dist (x3 , y3 ,  xl 9 ,  yl 9 )  ; 
=  node_dist (x3, y3, x20, y20) ; 
=  node_dist (x3 , y3 , x2 1 ,  y2 1 )  ; 

dl3; 
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d42  =  d24; 
d43  =  d34; 

d45  =  node_dist (x4 , y4 , x5 , y5 ) ; 
d46  =  node_dist (x4 , y4 , x6 , y6 ) ; 
d47  =  node_dist (x4 , y4  ,  x7  ,  y7  )  ; 
d4  8  =  node_dist (x4 ,  y4 ,  x8 ,  y8 )  ; 
d49  =  node_dist (x4 , y4 ,  x9 ,  y9 )  ; 
d410  =  node_dist (x4 , y4 ,  xlO ,  ylO )  ; 
d411  =  node_dist (x4 , y4  ,  xl 1 ,  yl 1 )  ; 
d412  =  node_dist (x4 , y4 , xl2  ,  yl2 )  ; 
d413  =  node_dist (x4, y4, xl3, yl3) ; 
d414  =  node_dist (x4  ,  y4  ,  xl4  ,  yl4  )  ; 
d415  =  node_dist (x4, y4, xl5, yl5) ; 
d416  =  node_dist (x4 ,  y4 ,  xl 6 ,  yl 6 )  ; 
d417  =  node_dist (x4 , y4 , xl7  ,  yl7 )  ; 
d418  =  node_dist (x4 , y4 ,  xl 8 ,  yl 8 )  ; 
d419  =  node_dist (x4 ,  y4 ,  xl 9 ,  yl 9 )  ; 
d420  =  node_dist (x4 , y4 ,  x20 ,  y20 )  ; 
d421  =  node_dist (x4 , y4 , x2 1 , y2 1 ) ; 

d51  =  dl5; 
d52  =  d25; 
d53  =  d35; 
d54  =  d45; 

d56  =  node_dist (x5 , y5 , x6 , y6 ) ; 
d57  =  node_dist (x5 , y5 ,  x7 ,  y7 )  ; 
d58  =  node_dist (x5 ,  y5 ,  x8 ,  y8 )  ; 
d59  =  node_dist (x5 , y5 ,  x9 ,  y9 )  ; 
dSlO  =  node_dist (x5, y5, xlO, ylO) ; 
dSll  =  node_dist (x5 , y5 , xl 1 ,  yl 1 ) ; 
d512  =  node_dist (x5 , y5 ,  xl2 ,  yl2 )  ; 
d513  =  node_dist (x5, y5, xl3, yl3) ; 
d514  =  node_dist (x5, y5, xl4, yl4) ; 
d515  =  node_dist (x5, y5, xl5,  yl5) ; 
d516  =  node_dist (x5 , y5 ,  xl 6 ,  yl 6 )  ; 
d517  =  node_dist (x5, y5, xl7, yl7) ; 
d518  =  node_dist (x5 , y5 , xl 8 ,  yl 8 ) ; 
d519  =  node_dist (x5 , y5 ,  xl 9 ,  yl 9 )  ; 
d520  =  node_dist (x5, y5, x20, y20) ; 
d521  =  node_dist (x5 , y5 , x2 1 ,  y2 1 ) ; 

d61  =  dl6; 
d62  =  d26; 
d63  =  d36; 
d64  =  d46; 
d65  =  d56; 

d67  =  node_dist (x6 , y6 , x7 , y7 ) ; 
d68  =  node_dist (x6 ,  y6 ,  x8 ,  y8 )  ; 
d69  =  node_dist (x6 , y6 ,  x9 ,  y9 )  ; 
d610  =  node_dist (x6, y6, xlO, ylO) ; 
d611  =  node_dist (x6 , y6 , xl 1 ,  yl 1 ) ; 
d612  =  node_dist (x6 , y6 ,  xl2 ,  yl2 )  ; 
d613  =  node_dist (x6, y6, xl3, yl3) ; 
d614  =  node_dist (x6, y6, xl4, yl4) ; 
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d615  =  node_dist (x6, y6, xl5, yl5) ; 
d616  =  node_dist (x6 ,  y6 ,  xl 6 ,  yl 6 )  ; 
d617  =  node_dist (x6, y6, xl7, yl7) ; 
d618  =  node_dist (x6 , y6 , xl 8 ,  yl 8 ) ; 
d619  =  node_dist (x6 , y6 ,  xl 9 ,  yl 9 )  ; 
d620  =  node_dist (x6, y6, x20, y20) ; 
d621  =  node_dist (x6 , y6 , x2 1 ,  y2 1 ) ; 

d71  =  dl7; 
d72  =  6.21; 

613  =  631; 
d74  =  d47; 
d75  =  d57; 
d76  =  d67; 

d78  =  node_dist (x7 , y7 , x8 , y8 ) ; 
d79  =  node_dist (x7 , y7 , x9 , y9 ) ; 
d710  =  node_dist (x7 , y7 ,  xlO ,  ylO )  ; 
d711  =  node_dist (x7 , y7  ,  xl 1 ,  yl 1 )  ; 
d712  =  node_dist (x7 , y7 , xl2  ,  yl2 )  ; 
d713  =  node_dist (x7, y7, xl3, yl3) ; 
d714  =  node_dist (x7  ,  y7  ,  xl4  ,  yl4  )  ; 
d715  =  node_dist (x7, y7, xl5, yl5) ; 
d716  =  node_dist (x7 ,  y7 ,  xl 6 ,  yl 6 )  ; 
d717  =  node_dist (x7 , y7 , xl7  ,  yl7 )  ; 
d718  =  node_dist (x7 , y7 ,  xl 8 ,  yl 8 )  ; 
d719  =  node_dist (x7 ,  y7 ,  xl 9 ,  yl 9 )  ; 
d720  =  node_dist (x7 , y7 ,  x20 ,  y20 )  ; 
d721  =  node_dist (x7 , y7 , x2 1 , y2 1 ) ; 

d81  =  dl8; 

d82  =  d28; 

d83  =  d38; 

d84  =  d48; 

d85  =  d58; 

d86  =  d68; 

d87  =  d78; 

d89  =  node_dist (x8 , y8 , x9 , y9 ) ; 

d810  =  node_dist (x8 , y8 , xlO , ylO ) ; 
d811  =  node_dist (x8 , y8 ,  xl 1 ,  yl 1 )  ; 
d812  =  node_dist (x8 , y8 ,  xl2 ,  yl2 )  ; 
d813  =  node_dist (x8, y8, xl3, yl3) ; 
d814  =  node_dist (x8 ,  y8 ,  xl4 ,  yl4 )  ; 
d815  =  node_dist (x8, y8, xl5, yl5) ; 
d816  =  node_dist (x8 ,  y8 ,  xl 6 ,  yl 6 )  ; 
d817  =  node_dist (x8 , y8 ,  xl7 ,  yl7 )  ; 
d818  =  node_dist (x8 , y8 ,  xl 8 ,  yl 8 )  ; 
d819  =  node_dist (x8 ,  y8 ,  xl 9 ,  yl 9 )  ; 
d820  =  node_dist (x8 , y8 ,  x20 ,  y20 )  ; 
d821  =  node_dist (x8 ,  y8 ,  x2 1 ,  y2 1 )  ; 

d91  =  dl9; 
d92  =  d29; 
d93  =  d39; 
d94  =  d49; 
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d95  =  d59; 
d96  =  d69; 
d97  =  d79; 
d98  =  d89; 

d910  =  node_dist (x9, y9, xlO, ylO) ; 
d911  =  node_dist (x9 ,  y9 ,  xl 1 ,  yl 1 )  ; 
d912  =  node_dist (x9 , y9 ,  xl2 ,  yl2 )  ; 
d913  =  node_dist (x9, y9, xl3, yl3) ; 
d914  =  node_dist (x9, y9, xl4, yl4) ; 
d915  =  node_dist (x9, y9, xl5,  yl5) ; 
d916  =  node_dist (x9 , y9 ,  xl 6 ,  yl 6 )  ; 
d917  =  node_dist (x9, y9, xl7, yl7) ; 
d918  =  node_dist (x9 , y9 , xl 8 ,  yl 8 ) ; 
d919  =  node_dist (x9 , y9 ,  xl 9 ,  yl 9 )  ; 
d920  =  node_dist (x9, y9, x20, y20) ; 
d921  =  node_dist (x9 , y9 , x2 1 ,  y2 1 ) ; 

dlOl  =  dllO; 
dl02  =  d210; 
dl03  =  d310; 
dl04  =  d410; 
dl05  =  dSlO; 
dl06  =  d610; 
dl07  =  d710; 
dl08  =  d810; 
dl09  =  d910; 

dlOll  =  node_dist (xlO, ylO, xll, yll) ; 
dl012  =  node_dist (xlO ,  ylO ,  xl2 ,  yl2 )  ; 
dl013  =  node_dist (xlO, ylO, xl3, yl3) ; 
dl014  =  node_dist (xlO , ylO  ,  xl4  ,  yl4  )  ; 
dlOlS  =  node_dist (xlO, ylO, xl5, yl5) ; 
dlOlG  =  node_dist (xlO, ylO, xl6, yl6) ; 
dl017  =  node_dist (xlO , ylO ,  xl7 ,  yl7 )  ; 
dl018  =  node_dist (xlO, ylO, xl8, yl8) ; 
dl019  =  node_dist (xlO, ylO, xl9, yl9) ; 
dl020  =  node_dist (xlO , ylO ,  x20 ,  y20 )  ; 
dl021  =  node_dist (xlO, ylO, x21, y21) ; 

dill  =  dill; 
dll2  =  d211; 
dll3  =  d311; 
dll4  =  d411; 
dll5  =  dSll; 
dll6  =  d611; 
dll7  =  d711; 
dll  8  =  d811; 
dll9  =  d911; 
dlllO  =  dlOll; 

dlll2  =  node_dist (xl 1 , yl 1 , xl2 , yl2 ) ; 
dlll3  =  node_dist (xll, yll, xl3,  yl3)  ; 
dlll4  =  node_dist (xll, yll, xl4, yl4) ; 
dlllS  =  node_dist (xll,  yll,  xl5,  yl5)  ; 
dlllG  =  node_dist (xl 1 , yl 1 ,  xl 6 ,  yl 6 )  ; 
dlll7  =  node_dist (xll, yll, xl7, yl7) ; 
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dlll8  =  node_dist (xl 1 ,  yl 1 ,  xl 8 ,  yl 8 )  ; 
dlll9  =  node_dist (xl 1 , yl 1 ,  xl 9 ,  yl 9 )  ; 
dll20  =  node_dist (xll, yll, x20, y20) ; 
dll21  =  node_dist (xl 1 , yl 1 ,  x2 1 ,  y2 1 )  ; 

dl21  =  dll2; 
dl22  =  d212; 
dl23  =  d312; 
dl24  =  d412; 
dl25  =  d512; 
dl26  =  d612; 
dl27  =  d712; 
dl28  =  d812; 
dl29  =  d912; 
dl210  =  dl012; 
dl211  =  dlll2; 

dl213  =  node_dist (xl2, yl2, xl3, yl3) ; 
dl214  =  node_dist (xl2 , yl2  ,  xl4  ,  yl4 )  ; 
dl215  =  node_dist (xl2, yl2, xl5, yl5) ; 
dl216  =  node_dist (xl2 ,  yl2 ,  xl 6 ,  yl 6 )  ; 
dl217  =  node_dist (xl2 , yl2  ,  xl7  ,  yl7  )  ; 
dl218  =  node_dist (xl2 , yl2 ,  xl 8 ,  yl 8 )  ; 
dl219  =  node_dist (xl2 , yl2 ,  xl 9 ,  yl 9 )  ; 
dl220  =  node_dist (xl2 , yl2 ,  x20 ,  y20 )  ; 
dl221  =  node_dist (xl2  ,  yl2  ,  x2 1 ,  y2 1 )  ; 

dl31  =  dll3; 
dl32  =  d213; 
dl33  =  d313; 
dl34  =  d413; 
dl35  =  d513; 
dl36  =  d613; 
dl37  =  d713; 
dl38  =  d813; 
dl39  =  d913; 
dl310  =  dl013; 
dl311  =  dlll3; 
dl312  =  dl213; 

dl314  =  node_dist (xl3, yl3, xl4, yl4) ; 
dl315  =  node_dist (xl3, yl3,  xl5,  yl5)  ; 
dl316  =  node_dist (xl3, yl3, xl6, yl6) ; 
dl317  =  node_dist (xl3, yl3, xl7, yl7) ; 
dl318  =  node_dist (xl3, yl3,  xl8,  yl8)  ; 
dl319  =  node_dist (xl3, yl3, xl9, yl9) ; 
dl320  =  node_dist (xl3, yl3, x20, y20) ; 
dl321  =  node_dist (xl3, yl3,  x21,  y21)  ; 

dl41  =  dll4; 
dl42  =  d214; 
dl43  =  d314; 
dl44  =  d414; 
dl45  =  d514; 
dl46  =  d614; 
dl47  =  d714; 
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dl48  = 

dl49  = 

dl410 

dl411 

dl412 

dl413 

dl415 

dl416 

dl417 

dl418 

dl419 

dl420 

dl421 

dl51  = 

dl52  = 

dl53  = 

dl54  = 

dl56  = 

dl57  = 

dl58  = 

dl59  = 

dlSlO 

dlSll 

dl512 

dl513 

dl514 

dl516 

dl517 

dl518 

dl519 

dl520 

dl521 

dl61  = 
dl62  = 
dl63  = 
dl64  = 

dies  = 

diee  = 

dl67  = 

dl68  = 

dl69  = 

dieio 

dieil 

dl612 

dl613 

dl614 

dieis 

dl617 

dl618 

dl619 

dl620 

dl621 


d814; 

d914; 

=  dl014; 

=  dlll4; 

=  dl214; 

=  dl314; 

=  node_dist (xl4, yl4, xl5, yl5) ; 
=  node_dist (xl4, yl4, xie, yl6) ; 
=  node_dist (xl4 , yl4  ,  xl7  ,  yl7  )  ; 
=  node_dist (xl4, yl4, xl8, yl8) ; 
=  node_dist (xl4, yl4, xl9, yl9) ; 
=  node_dist (xl4 , yl4 ,  x20 ,  y20 )  ; 
=  node_dist (xl4, yl4, x21, y21) ; 

dllS; 

d215; 

d315; 

d415; 

d615; 

d715; 

d815; 

d915; 

=  dlOlS; 

=  dlllS; 

=  dl215; 

=  dl315; 

=  dl415; 

=  node_dist (xl5, yl5, xie,  yl6)  ; 
=  node_dist (xl5, yl5, xl7, yl7) ; 
=  node_dist (xl5, yl5, xl8, yl8) ; 
=  node_dist (xl5, yl5, xl9, yl9) ; 
=  node_dist (xl5, yl5,  x20,  y20)  ; 
=  node_dist (xl5, yl5, x21, y21) ; 

dll6; 
d2 1 6  ; 
d316; 
d4 1 6  ; 
d516; 
d616; 
d7 1 6  ; 
d816; 
d916; 

=  dl016; 

=  dlll6; 

=  dl216; 

=  dl316; 

=  dl416; 

=  dl516; 

=  node_dist (xie, yie, xl7,  yl7)  ; 
=  node_dist (xl e , yl e , xl 8 , yl 8 ) ; 
=  node_dist (xl e , yl e , xl 9 ,  yl 9 ) ; 
=  node_dist (xie, yie, x20, y20) ; 
=  node_dist (xie, yie,  x20,  y21)  ; 
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dl71  = 

dll7; 

dl72  = 

d217; 

dl73  = 

d317; 

dl74  = 

d417; 

dl75  = 

d517; 

dl76  = 

d617; 

dl77  = 

d717; 

dl78  = 

d817; 

dl79  = 

d917; 

dl710 

=  dl017; 

dl711 

=  dlll7; 

dl712 

=  dl217; 

dl713 

=  dl317; 

dl714 

=  dl417; 

dl715 

=  dl517; 

dl716 

=  dl617; 

dl718 

=  node  dist (xl7, yl7, xl8, yl8) 

dl719 

=  node  dist (xl7, yl7, xl9, yl9) 

dl720 

=  node  dist (xl7  ,  yl7  ,  x20  ,  y20  ) 

dl721 

=  node  dist (xl7, yl7, x21, y21) 

dl81  = 

dll8; 

dl82  = 

d2 1 8  ; 

dl83  = 

d318; 

dl84  = 

d4 1 8  ; 

dl85  = 

d518; 

dl8  6  = 

d618; 

dl87  = 

d7 1 8  ; 

dl8  8  = 

d818; 

dl8  9  = 

d918; 

dl810 

=  dl018; 

dl811 

=  dlll8; 

dl812 

=  dl218; 

dl813 

=  dl318; 

dl814 

=  dl418; 

dl815 

=  dl518; 

dl816 

=  dl618; 

dl817 

=  dl718; 

dl819 

=  node  dist (xl 8 , yl 8 , xl 9 , yl 9 ) 

dl820 

=  node  dist (xl8, yl8, x20, y20) 

dl821 

=  node  dist (xl 8 , yl 8 , x2 1 , y2 1 ) 

dl91  = 

dll9; 

dl92  = 

d2 1 9  ; 

dl93  = 

d319; 

dl94  = 

d4 1 9  ; 

dl95  = 

d519; 

dl96  = 

d619; 

dl97  = 

d7 1 9  ; 

dl98  = 

d819; 

dl99  = 

d919; 

dl910 

=  dl019; 

dl911 

=  dlll9; 

dl912 

=  dl219; 

dl913 

=  dl319; 

dl914 

=  dl419; 

dl915 

=  dl519; 

dl916 

=  dl619; 

dl917 

=  dl719; 

dl918 

=  dl819; 

dl920 

=  node  dist 

dl921 

=  node  dist 

d201  = 

dl20; 

d202  = 

d220; 

d203  = 

d32  0; 

d204  = 

d420; 

d205  = 

d52  0; 

d206  = 

d62  0; 

d207  = 

d720; 

d208  = 

d820; 

d209  = 

d92  0; 

d2010 

=  dl020; 

d2011 

=  dll20; 

d2012 

=  dl220; 

d2013 

=  dl320; 

d2014 

=  dl420; 

d2015 

=  dl520; 

d2016 

=  dl620; 

d2017 

=  dl720; 

d2018 

=  dl820; 

d2019 

=  dl920; 

d2021 

=  node  dist 

d211  = 

dl21; 

d212  = 

d221; 

d213  = 

d321; 

d214  = 

d421; 

d215  = 

d521; 

d216  = 

d621; 

d217  = 

d721; 

d218  = 

d821; 

d219  = 

d921; 

d2110 

=  dl021; 

d2111 

=  dll21; 

d2112 

=  dl221; 

d2113 

=  dl321; 

d2114 

=  dl421; 

d2115 

=  dl521; 

d2116 

=  dl621; 

d2117 

=  dl721; 

d2118 

=  dl821; 

d2119 

=  dl921; 

d2120 

=  d2021; 

%Calculate  the  recieved  power  due  to  freespace 


losses 
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Prxl2  =  pathlosscalc ( Ptx, dl2 ) ; 
Prxl3  =  pathlosscalc (Ptx, dl3) ; 
Prxl4  =  pathlosscalc (Ptx, dl4) ; 
PrxlS  =  pathlosscalc (Ptx, dl5) ; 
Prxl6  =  pathlosscalc ( Ptx, dl 6 ) ; 
Prxl7  =  pathlosscalc (Ptx, dl7) ; 
PrxlS  =  pathlosscalc ( Ptx, dl 8 ) ; 
Prxl9  =  pathlosscalc ( Ptx, dl 9 ) ; 
PrxllO  =  pathlosscalc (Ptx, dllO)  ; 
Prxlll  =  pathlosscalc ( Ptx, dl 1 1 ) ; 
Prxll2  =  pathlosscalc ( Ptx, dl 12 ) ; 
Prxll3  =  pathlosscalc (Ptx, dll3) ; 
Prxll4  =  pathlosscalc (Ptx, dll4) ; 
PrxllS  =  pathlosscalc (Ptx, dll5) ; 
PrxllG  =  pathlosscalc ( Ptx, dl 1 6 ) ; 
Prxll7  =  pathlosscalc (Ptx, dll7)  ; 
PrxllS  =  pathlosscalc ( Ptx, dl 1 8 ) ; 
Prxll9  =  pathlosscalc ( Ptx, dl 1 9 ) ; 
Prxl20  =  pathlosscalc (Ptx, dl20)  ; 
Prxl21  =  pathlosscalc ( Ptx, dl2 1 ) ; 

Prx21  =  pathlosscalc ( Ptx, d2 1 ) ; 
Prx23  =  pathlosscalc (Ptx, d23) ; 
Prx24  =  pathlosscalc (Ptx, d24) ; 
Prx25  =  pathlosscalc (Ptx, d25) ; 
Prx26  =  pathlosscalc ( Ptx, d2 6 ) ; 
Prx27  =  pathlosscalc (Ptx, d27) ; 
Prx28  =  pathlosscalc ( Ptx, d2 8 )  ; 
Prx29  =  pathlosscalc ( Ptx, d2 9 ) ; 
Prx210  =  pathlosscalc (Ptx, d210)  ; 
Prx211  =  pathlosscalc ( Ptx, d2 1 1 ) ; 
Prx212  =  pathlosscalc ( Ptx, d2 12 ) ; 
Prx213  =  pathlosscalc (Ptx, d213) ; 
Prx214  =  pathlosscalc (Ptx, d214) ; 
Prx215  =  pathlosscalc (Ptx, d215) ; 
Prx216  =  pathlosscalc ( Ptx, d2 1 6 ) ; 
Prx217  =  pathlosscalc (Ptx, d217) ; 
Prx218  =  pathlosscalc ( Ptx,  d2 1 8 )  ; 
Prx219  =  pathlosscalc ( Ptx, d2 1 9 ) ; 
Prx220  =  pathlosscalc (Ptx, d220) ; 
Prx221  =  pathlosscalc ( Ptx, d22 1 )  ; 

Prx31  =  pathlosscalc (Ptx, d31) ; 
Prx32  =  pathlosscalc ( Ptx, d32 ) ; 
Prx34  =  pathlosscalc ( Ptx, d34 ) ; 
Prx35  =  pathlosscalc ( Ptx, d35 ) ; 
Prx36  =  pathlosscalc (Ptx, d36) ; 
Prx37  =  pathlosscalc ( Ptx, d37 ) ; 
Prx38  =  pathlosscalc (Ptx, d38) ; 
Prx39  =  pathlosscalc (Ptx, d39) ; 
Prx310  =  pathlosscalc (Ptx,  d310)  ; 
Prx311  =  pathlosscalc (Ptx, d311) ; 
Prx312  =  pathlosscalc (Ptx, d312) ; 
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Prx313 
Prx314 
Prx315 
Prx316 
Prx317 
Prx318 
Prx319 
Prx320 
Prx32 1 


pathlosscalc (Ptx, d313) 
pathlosscalc (Ptx, d314) 
pathlosscalc (Ptx, d315) 
pathlosscalc (Ptx, d316) 
pathlosscalc (Ptx, d317) 
pathlosscalc (Ptx, d318) 
pathlosscalc (Ptx, d319) 
pathlosscalc (Ptx, d320) 
pathlosscalc ( Ptx, d32 1 ) 


Prx4 1 
Prx42 
Prx43 
Prx45 
Prx4  6 
Prx47 
Prx4  8 
Prx4  9 
Prx410 
Prx4 1 1 
Prx4 12 
Prx413 
Prx414 
Prx415 
Prx4 1 6 
Prx417 
Prx4 1 8 
Prx4 1 9 
Prx420 
Prx42 1 


pathlosscalc ( Ptx,  d4 1 )  ; 
pathlosscalc ( Ptx,  d42  )  ; 
pathlosscalc (Ptx,  d43)  ; 
pathlosscalc (Ptx,  d45)  ; 
pathlosscalc ( Ptx, d4 6 )  ; 
pathlosscalc (Ptx,  d4  7)  ; 
pathlosscalc ( Ptx,  d4 8 )  ; 
pathlosscalc ( Ptx,  d4 9 )  ; 

^  pathlosscalc (Ptx, d410) 
^  pathlosscalc ( Ptx, d4 1 1 ) 
^  pathlosscalc ( Ptx, d4 12 ) 
^  pathlosscalc (Ptx, d413) 
^  pathlosscalc (Ptx, d414) 
^  pathlosscalc (Ptx, d415) 
^  pathlosscalc ( Ptx, d4 1 6 ) 
^  pathlosscalc (Ptx, d417) 
^  pathlosscalc ( Ptx, d4 1 8 ) 
^  pathlosscalc ( Ptx, d4 1 9 ) 
^  pathlosscalc (Ptx, d420) 
^  pathlosscalc ( Ptx, d42 1 ) 


PrxSl 

Prx52 

Prx53 

Prx54 

Prx56 

Prx57 

Prx58 

Prx59 

PrxSlO 

PrxSll 

Prx512 

Prx513 

Prx514 

Prx515 

Prx516 

Prx517 

Prx518 

Prx519 

Prx520 

Prx52 1 


pathlosscalc (Ptx,  dSl)  ; 
pathlosscalc ( Ptx,  d52  )  ; 
pathlosscalc ( Ptx,  d53 )  ; 
pathlosscalc ( Ptx,  d54  )  ; 
pathlosscalc (Ptx, d56)  ; 
pathlosscalc ( Ptx,  d57  )  ; 
pathlosscalc (Ptx,  d58)  ; 
pathlosscalc (Ptx, d59)  ; 

^  pathlosscalc (Ptx, dSlO) 
^  pathlosscalc (Ptx,  dSll) 
^  pathlosscalc (Ptx,  d512) 
^  pathlosscalc (Ptx, d513) 
^  pathlosscalc (Ptx,  d514) 
^  pathlosscalc (Ptx,  d515) 
^  pathlosscalc (Ptx, d516) 
^  pathlosscalc (Ptx,  d517) 
^  pathlosscalc (Ptx,  d518) 
^  pathlosscalc (Ptx, d519) 
^  pathlosscalc (Ptx,  d520) 
^  pathlosscalc ( Ptx,  d52 1 ) 


PrxGl  =  PrxlG; 
Prx62  =  Prx26; 


Prx63  =  Prx36; 

Prx64  =  Prx46; 

Prx65  =  Prx56; 

Prx67  =  pathlosscalc ( Ptx, d67 ) ; 
Prx68  =  pathlosscalc (Ptx,  d68)  ; 
Prx69  =  pathlosscalc (Ptx, d69) ; 
PrxGlO  =  pathlosscalc (Ptx, d610) ; 
PrxGll  =  pathlosscalc (Ptx, d611)  ; 
Prx612  =  pathlosscalc (Ptx, d612) ; 
Prx613  =  pathlosscalc (Ptx, d613) ; 
Prx614  =  pathlosscalc (Ptx, d614)  ; 
Prx615  =  pathlosscalc (Ptx, d615) ; 
Prx616  =  pathlosscalc (Ptx, d616) ; 
Prx617  =  pathlosscalc (Ptx, d617)  ; 
Prx618  =  pathlosscalc (Ptx, d618) ; 
Prx619  =  pathlosscalc (Ptx, d619) ; 
Prx620  =  pathlosscalc (Ptx, d620)  ; 
Prx621  =  pathlosscalc ( Ptx, d62 1 ) ; 

Prx71  =  Prxl7; 

Prx72  =  Prx27; 

Prx73  =  Prx37; 

Prx74  =  Prx47; 

Prx75  =  Prx57; 

Prx76  =  Prx67; 

Prx78  =  pathlosscalc ( Ptx, d7 8 ) ; 
Prx79  =  pathlosscalc ( Ptx, d7 9 ) ; 
Prx710  =  pathlosscalc (Ptx, d710) ; 
Prx711  =  pathlosscalc ( Ptx, d7 1 1 ) ; 
Prx712  =  pathlosscalc ( Ptx, d7 12 ) ; 
Prx713  =  pathlosscalc (Ptx, d713) ; 
Prx714  =  pathlosscalc (Ptx, d714)  ; 
Prx715  =  pathlosscalc (Ptx, d715) ; 
Prx716  =  pathlosscalc ( Ptx, d7 1 6 ) ; 
Prx717  =  pathlosscalc (Ptx, d717)  ; 
Prx718  =  pathlosscalc ( Ptx, d7 1 8 ) ; 
Prx719  =  pathlosscalc ( Ptx, d7 1 9 ) ; 
Prx720  =  pathlosscalc (Ptx, d720)  ; 
Prx721  =  pathlosscalc (Ptx, d721) ; 

Prx81  =  Prxl8; 

Prx82  =  Prx28; 

Prx83  =  Prx38; 

Prx84  =  Prx48; 

Prx85  =  Prx58; 

Prx86  =  Prx68; 

Prx87  =  Prx78; 

Prx89  =  pathlosscalc ( Ptx, d8 9 ) ; 
Prx810  =  pathlosscalc (Ptx, d810) ; 
Prx811  =  pathlosscalc ( Ptx, d8 1 1 ) ; 
Prx812  =  pathlosscalc ( Ptx, d8 12 ) ; 
Prx813  =  pathlosscalc (Ptx, d813) ; 
Prx814  =  pathlosscalc (Ptx, d814) ; 
Prx815  =  pathlosscalc (Ptx, d815) ; 
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Prx816  =  pathlosscalc ( Ptx, d8 1 6 ) ; 
Prx817  =  pathlosscalc (Ptx, d817)  ; 
Prx818  =  pathlosscalc ( Ptx, d8 1 8 ) ; 
Prx819  =  pathlosscalc ( Ptx, d8 1 9 ) ; 
Prx820  =  pathlosscalc (Ptx, d820)  ; 
Prx821  =  pathlosscalc ( Ptx, d82 1 ) ; 


Prx91  = 
Prx92  = 
Prx93  = 
Prx94  = 
Prx95  = 
Prx96  = 
Prx97  = 
Prx98  = 
Prx910  = 
Prx911  = 
Prx912  = 
Prx913  = 
Prx914  = 
Prx915  = 
Prx916  = 
Prx917  = 
Prx918  = 
Prx919  = 
Prx920  = 
Prx921  = 


Prxl  9 ; 

Prx2  9 ; 

Prx39; 

Prx4  9 ; 

Prx59; 

Prx69; 

Prx7  9 ; 

Prx8  9 ; 

pathlosscalc (Ptx,  d910)  ; 
pathlosscalc (Ptx,  d911)  ; 
pathlosscalc (Ptx,  d912)  ; 
pathlosscalc (Ptx,  d913)  ; 
pathlosscalc (Ptx, d914)  ; 
pathlosscalc (Ptx,  d915)  ; 
pathlosscalc (Ptx,  d916)  ; 
pathlosscalc (Ptx,  d917)  ; 
pathlosscalc (Ptx,  d918)  ; 
pathlosscalc (Ptx,  d919)  ; 
pathlosscalc (Ptx,  d92  0)  ; 
pathlosscalc ( Ptx,  d92 1 )  ; 


PrxlOl  = 
Prxl02  = 
Prxl03  = 
Prxl04  = 
PrxlOS  = 
PrxlOG  = 
Prxl07  = 
Prxl08  = 
Prxl09  = 
PrxlOll  = 
Prxl012  = 
Prxl013  = 
Prxl014  = 
PrxlOlS  = 
PrxlOlG  = 
Prxl017  = 
Prxl018  = 
Prxl019  = 
Prxl020  = 
Prxl021  = 


Prxl  10; 

Prx210; 

Prx310 ; 

Prx410; 

PrxSlO ; 

PrxGlO ; 

Prx710; 

Prx810; 

Prx910 ; 

pathlosscalc (Ptx,  dlOll)  ; 
pathlosscalc (Ptx, dl012)  ; 
pathlosscalc (Ptx, dl013)  ; 
pathlosscalc (Ptx, dl014)  ; 
pathlosscalc (Ptx,  dlOlS)  ; 
pathlosscalc (Ptx, dlOlG)  ; 
pathlosscalc (Ptx,  dl017)  ; 
pathlosscalc (Ptx,  dl018)  ; 
pathlosscalc (Ptx, dl019)  ; 
pathlosscalc (Ptx,  dl02  0)  ; 
pathlosscalc (Ptx,  dl021)  ; 


Prxl 11  =  Prxl 11; 
Prxl 12  =  Prx211; 
Prxl 13  =  Prx311; 
Prxl 14  =  Prx411; 
Prxl 15  =  Prx511; 
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PrxllG  =  PrxGll; 

Prxll7  =  Prx711; 

PrxllS  =  PrxSll; 

Prxll9  =  Prx911; 

PrxlllO  =  PrxlOll; 

Prxlll2  =  pathlosscalc ( Ptx, dl 1 12 ) ; 
PrxlllS  =  pathlosscalc (Ptx,  dlll3)  ; 
Prxlll4  =  pathlosscalc (Ptx, dlll4) ; 
PrxlllS  =  pathlosscalc (Ptx, dlllS) ; 
PrxlllG  =  pathlosscalc ( Ptx, dl 1 1 6 )  ; 
Prxlll7  =  pathlosscalc (Ptx, dlll7) ; 
PrxlllS  =  pathlosscalc ( Ptx, dl 1 1 8 ) ; 
Prxlll9  =  pathlosscalc ( Ptx, dl 1 1 9 ) ; 
Prxll20  =  pathlosscalc (Ptx, dll20) ; 
Prxll21  =  pathlosscalc ( Ptx, dl 12 1 ) ; 

Prxl21  =  Prxll2; 

Prxl22  =  Prx212; 

Prxl23  =  Prx312; 

Prxl24  =  Prx412; 

Prxl25  =  Prx512; 

Prxl26  =  Prx612; 

Prxl27  =  Prx712; 

Prxl28  =  Prx812; 

Prxl29  =  Prx912; 

Prxl210  =  Prxl012; 

Prxl211  =  Prxlll2; 

Prxl213  =  pathlosscalc (Ptx, dl213) ; 
Prxl214  =  pathlosscalc (Ptx, dl214) ; 
Prxl215  =  pathlosscalc (Ptx, dl215) ; 
Prxl216  =  pathlosscalc ( Ptx,  dl2 1 6 )  ; 
Prxl217  =  pathlosscalc (Ptx, dl217) ; 
Prxl218  =  pathlosscalc ( Ptx, dl2 1 8 ) ; 
Prxl219  =  pathlosscalc ( Ptx, dl2 1 9 )  ; 
Prxl220  =  pathlosscalc (Ptx, dl220) ; 
Prxl221  =  pathlosscalc ( Ptx, dl22 1 ) ; 

Prxl31  =  Prxll3; 

Prxl32  =  Prx213; 

Prxl33  =  Prx313; 

Prxl34  =  Prx413; 

Prxl35  =  Prx513; 

Prxl36  =  Prx613; 

Prxl37  =  Prx713; 

Prxl38  =  Prx813; 

Prxl39  =  Prx913; 

Prxl310  =  Prxl013; 

Prxl311  =  Prxlll3; 

Prxl312  =  Prxl213; 

Prxl314  =  pathlosscalc (Ptx, dl314) ; 
Prxl315  =  pathlosscalc (Ptx, dl315) ; 
Prxl316  =  pathlosscalc (Ptx, dl316) ; 
Prxl317  =  pathlosscalc (Ptx,  dl317)  ; 
Prxl318  =  pathlosscalc (Ptx, dl318) ; 
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Prxl319  =  pathlosscalc (Ptx, dl319) ; 
Prxl320  =  pathlosscalc (Ptx, dl320) ; 
Prxl321  =  pathlosscalc (Ptx, dl321)  ; 


Prxl41  = 
Prxl42  = 
Prxl43  = 
Prxl44  = 
Prxl45  = 
Prxl46  = 
Prxl47  = 
Prxl48  = 
Prxl49  = 
Prxl410  = 
Prxl411  = 
Prxl412  = 
Prxl413  = 
Prxl415  = 
Prxl416  = 
Prxl417  = 
Prxl418  = 
Prxl419  = 
Prxl420  = 
Prxl421  = 


Prxll4; 

Prx214; 

Prx314 ; 

Prx414; 

Prx514  ; 

Prx614 ; 

Prx714; 

Prx814; 

Prx914 ; 

Prxl014; 

Prxlll4; 

Prxl214; 

Prxl314 ; 

pathlosscalc (Ptx, dl415) ; 
pathlosscalc (Ptx, dl416)  ; 
pathlosscalc (Ptx,  dl417)  ; 
pathlosscalc (Ptx,  dl418)  ; 
pathlosscalc (Ptx, dl419)  ; 
pathlosscalc (Ptx,  dl42  0)  ; 
pathlosscalc (Ptx,  dl421)  ; 


PrxlSl  = 
Prxl52  = 
Prxl53  = 
Prxl54  = 
Prxl55  = 
Prxl56  = 
Prxl57  = 
Prxl58  = 
Prxl59  = 
PrxlSlO  = 
PrxlSll  = 
Prxl512  = 
Prxl513  = 
Prxl514  = 
PrxlSlG  = 
Prxl517  = 
Prxl518  = 
Prxl519  = 
Prxl520  = 
Prxl521  = 


PrxllS; 

Prx215; 

Prx315; 

Prx415; 

Prx515; 

Prx615; 

Prx715; 

Prx815; 

Prx915; 

PrxlOlS; 

PrxlllS; 

Prxl215; 

Prxl315; 

Prxl415; 

pathlosscalc (Ptx, dl516)  ; 
pathlosscalc (Ptx, dl517) ; 
pathlosscalc (Ptx, dl518) ; 
pathlosscalc (Ptx, dl519) ; 
pathlosscalc (Ptx, dl52  0)  ; 
pathlosscalc (Ptx, dl521)  ; 


PrxlGl  =  PrxllG; 
Prxl62  =  Prx216; 
Prxl63  =  Prx316; 
Prxl64  =  Prx416; 
Prxl65  =  Prx516; 
Prxl66  =  Prx616; 
Prxl67  =  Prx716; 
Prxl68  =  Prx816; 
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Prxl69  =  Prx916; 

PrxlGlO  =  PrxlOlG; 

PrxlGll  =  PrxlllG; 

Prxl612  =  Prxl216; 

Prxl613  =  Prxl316; 

Prxl614  =  Prxl416; 

Prxl615  =  Prxl516; 

Prxl617  =  pathlosscalc (Ptx, dl617) ; 
Prxl618  =  pathlosscalc (Ptx, dl618) ; 
Prxl619  =  pathlosscalc (Ptx, dl619) ; 
Prxl620  =  pathlosscalc (Ptx, dl620) ; 
Prxl621  =  pathlosscalc ( Ptx, dl 62 1 )  ; 

Prxl71  =  Prxll7; 

Prxl72  =  Prx217; 

Prxl73  =  Prx317; 

Prxl74  =  Prx417; 

Prxl75  =  Prx517; 

Prxl76  =  Prx617; 

Prxl77  =  Prx717; 

Prxl78  =  Prx817; 

Prxl79  =  Prx917; 

Prxl710  =  Prxl017; 

Prxl711  =  Prxlll7; 

Prxl712  =  Prxl217; 

Prxl713  =  Prxl317; 

Prxl714  =  Prxl417; 

Prxl715  =  Prxl517; 

Prxl716  =  Prxl617; 

Prxl718  =  pathlosscalc (Ptx, dl718) ; 
Prxl719  =  pathlosscalc (Ptx,  dl719)  ; 
Prxl720  =  pathlosscalc (Ptx, dl720) ; 
Prxl721  =  pathlosscalc (Ptx, dl721) ; 

Prxl81  =  Prxll8; 

Prxl82  =  Prx218; 

Prxl83  =  Prx318; 

Prxl84  =  Prx418; 

Prxl85  =  Prx518; 

Prxl86  =  Prx618; 

Prxl87  =  Prx718; 

Prxl88  =  Prx818; 

Prxl89  =  Prx918; 

Prxl810  =  Prxl018; 

Prxl811  =  Prxlll8; 

Prxl812  =  Prxl218; 

Prxl813  =  Prxl318; 

Prxl814  =  Prxl418; 

Prxl815  =  Prxl518; 

Prxl816  =  Prxl618; 

Prxl817  =  Prxl718; 

Prxl819  =  pathlosscalc ( Ptx, dl 8 1 9 ) ; 
Prxl820  =  pathlosscalc (Ptx, dl820) ; 
Prxl821  =  pathlosscalc ( Ptx, dl 82 1 ) ; 
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Prxl91  = 

Prxl 1 9 ; 

Prxl92  = 

Prx2 19; 

Prxl93  = 

Prx319; 

Prxl94  = 

Prx4 19; 

Prxl95  = 

Prx519; 

Prxl96  = 

Prx619; 

Prxl97  = 

Prx7 19; 

Prxl98  = 

Prx8 1 9 ; 

Prxl99  = 

Prx919; 

Prxl910 

=  Prxl019; 

Prxl911 

=  Prxlll9; 

Prxl912 

=  Prxl219; 

Prxl913 

=  Prxl319; 

Prxl914 

=  Prxl419; 

Prxl915 

=  Prxl519; 

Prxl916 

=  Prxl 61 9; 

Prxl 917 

=  Prxl719; 

Prxl918 

=  Prxl819; 

Prxl920 

=  pathlosscalc (Ptx, dl920)  ; 

Prxl 92 1 

=  pathlosscalc ( Ptx, dl 92 1 ) ; 

Prx201  = 

Prxl 20 ; 

Prx202  = 

Prx220 ; 

Prx203  = 

Prx320; 

Prx204  = 

Prx420 ; 

Prx205  = 

Prx52  0 ; 

Prx206  = 

Prx62  0 ; 

Prx207  = 

Prx720 ; 

Prx208  = 

Prx820 ; 

Prx209  = 

Prx920 ; 

Prx2010 

=  Prxl 020; 

Prx2011 

=  Prxll20; 

Prx2012 

=  Prxl220; 

Prx2013 

=  Prxl 320; 

Prx2014 

=  Prxl 420; 

Prx2015 

=  Prxl 520; 

Prx2016 

=  Prxl 620; 

Prx2017 

=  Prxl 720; 

Prx2018 

=  Prxl 820; 

Prx2019 

=  Prxl 920; 

Prx2021 

=  pathlosscalc (Ptx,  d2021)  ; 

Prx211  = 

Prxl 2 1 ; 

Prx212  = 

Prx22 1 ; 

Prx213  = 

Prx32 1 ; 

Prx214  = 

Prx42 1 ; 

Prx215  = 

Prx52 1 ; 

Prx216  = 

Prx62 1 ; 

Prx217  = 

Prx721; 

Prx218  = 

Prx82 1 ; 

Prx219  = 

Prx92 1 ; 

Prx2110 

=  Prxl 021; 

Prx2 111 

=  Prxl 121; 
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Prx2112  =  Prxl221 
Prx2113  =  Prxl321 
Prx2114  =  Prxl421 
Prx2115  =  Prxl521 
Prx2116  =  Prxl621 
Prx2117  =  Prxl721 
Prx2118  =  Prxl821 
Prx2119  =  Prxl921 
Prx2120  =  Prx2021 

%Calculate  Signal 

%SIR  for  transmis 

SIR2froml  =  Prxl2 
SIR6froml  =  Prxl6 
SIR7froml  =  Prxl7 

%SIR  for  transmis 

SIR3from2  =  Prx23 
SIR6from2  =  Prx23 
SIR7from2  =  Prx23 
SIR8from2  =  Prx23 

%SIR  for  transmis 

SIR7from3  =  Prx37 
SIR8from3  =  Prx38 
SIR9from3  =  Prx39 

%SIR  for  transmis 

SIR3from4  =  Prx43 
SIR5from4  =  Prx45 
SIR8from4  =  Prx48 
SIR9from4  =  Prx49 
SIR10from4  =  Prx4 

%SIR  for  transmis 

SIR4from5  =  Prx54 
SIR9from5  =  Prx59 
SIRlOfromS  =  Prx5 

%SIR  for  transmis 

SIR7from6  =  Prx67 
SIRllfromG  =  Prx6 
SIR12from6  =  Prx6 


to  Interference  Ratios 

sions  from  1  to  2,6,7 

/ (Prx32+Prx52+Prxll2+Prxl32+Prxl52) ; 

/ (Prx36tPrx56+Prxll6+Prxl36tPrxl56) ; 

/ (Prx37+Prx57+Prxll7+Prxl37+Prxl57) ; 

sions  from  2  to  3, 6, 7, 8 

/ (Prxl23+Prxl43tPrx43) ; 

/ (Prxl26+Prxl46tPrx46) ; 

/ (Prxl27+Prxl47+Prx47) ; 

/ (Prxl28+Prxl48+Prx48) ; 

sions  from  3  to  7,8,9 

/ (Prxl7+Prx57+Prxll7+Prxl37+Prxl57) ; 

/ (Prxl8+Prx58+Prxll8+Prxl38+Prxl58) ; 

/ (Prxl9tPrx59+Prxll9+Prxl39tPrxl59) ; 

sions  from  4  to  3,5,8,9,10 

/ (Prx23tPrxl23  +  Prxl43)  ; 

/  (Prx25tPrxl25  +  Prxl45)  ; 

/ (Prx28  +  Prxl28  +  Prxl48)  ; 

/ (Prx29+Prxl29+Prxl49) ; 

10/ (Prx210+Prxl210+Prxl410) ; 

sions  from  5  to  4,9,10 

/ (Prx34+Prxl4+Prxll4+Prxl34+Prxl54) ; 

/ (Prx39tPrxl9+Prxll9+Prxl39+Prxl59) ; 

10/ (Prx310+Prxll0tPrxlll0tPrxl310tPrxl510) 

sions  from  6  to  7,11,12 

/ (Prx87  +  Prxl07  +  Prxl67  +  Prxl87  +  Prx207)  ; 

1 1 /  ( Prx8 1 1  +  PrxlO 1 1  +  Prxl 61 1  +  Prxl 8 1 1  +  Prx2  0 1 1 
12 / ( Prx8 12+PrxlO 12+Prxl 612+Prxl 8 12+Prx20 12 
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SIR  for  transmissions  from  7  to  8,11,12,13 


SIR8from7  =  Prx87/ (Prx98  +  Prxl78  +  Prxl98)  ; 


SIRllfrom7  = 

Prx811 

SIR12from7  = 

Prx8 12 

SIR13from7  = 

Prx813 

%SIR  for  transmissii 

SIR12from8  = 

Prx8 12 

SIR13from8  = 

Prx813 

SIR14from8  = 

Prx814 

%SIR  for  transmissii 

SIR8from9  =  Prx98/ ( 

SIR13from9  = 

Prx913 

SIR14from9  = 

Prx914 

SIR15from9  = 

Prx915 

%SIR  for  tran 

smissii 

SIR9fromlO  = 

Prxl09 

SIR14fromlO  = 

PrxlO 

SIRlSfromlO  = 

PrxlO 

%SIR  for  tran 

smissii 

SIR12fromll  = 

Prxll 

SIRlGfromll  = 

Prxl  1 

SIR17fromll  = 

Prxll 

%SIR  for  tran 

smissii 

SIR13froml2  = 

Prxl2 

SIR16froml2  = 

Prxl2 

SIR17froml2  = 

Prxl2 

SIR18froml2  = 

Prxl2 

%SIR  for  tran 

smissii 

SIR17froml3  = 

Prxl3 

SIR18froml3  = 

Prxl3 

SIR19froml3  = 

Prxl3 

%SIR  for  tran 

smissii 

SIR13froml4  = 

Prxl4 

SIR18froml4  = 

Prxl4 

SIR19froml4  = 

Prxl4 

SIR20froml4  = 

Prxl4: 

to  12,13,14 


9  to  8, 13, 14, 15 


76 


%SIR  for  transmissions  from  15  to  14,19,20 

SIR14froml5  =  Prxl514/ (Prxll4+Prx314+Prx514+Prxlll4+Prxl314) ; 
SIR19froml5  =  Prxl519/ (Prxll9+Prx319+Prx519+Prxlll9+Prxl319) ; 
SIR20froml5  =  Prxl520/ (Prxl20  +  Prx320tPrx520tPrxll20tPrxl320)  ; 

%SIR  for  transmissions  from  16  to  17,21 

SIR17froml6  =  Prxl617/ (Prx617+Prx817+Prxl017+Prxl817+Prx2017) ; 
SIR21froml6  =  Prxl621/ (Prx621+Prx821+Prxl021+Prxl821+Prx2021) ; 


%SIR  for  transmissions  from  17  to  18,21 

SIR18froml7  =  Prxl718/ (Prx718  +  Prx918  +  Prxl918)  ; 
SIR21froml7  =  Prxl721/ (Prx721+Prx921+Prxl921) ; 


%SIR  for  transmissions  from 


18  to  21 


SIR21froml8  =  Prxl821/ (Prxl621  +  Prx2021  +  Prx621  +  Prx821  +  Prxl021)  ; 

%SIR  for  transmissions  from  19  to  18,21 

SIR18froml9  =  Prxl918/ (Prxl718+Prx718+Prx918) 

SIR21froml9  =  Prxl921/ (Prxl721+Prx721+Prx921) 


%SIR  for  transmissions  from  20  to  19,21 

SIR19from20  =  Prx2019/ (Prxl819+Prxl619tPrx619+Prx819tPrxl019) 
SIR21from20  =  Prx2021/ (Prxl821+Prxl621+Prx621+Prx821+Prxl021) 


W  =  ([SIR2froml  SIRGfroml  SIR7froml  SIR6from2  SIR7from2  SIR8from2  .. 
SIR3from2  SIR7from3  SIR8from3  SIR9from3  SIR8from4  SIR9from4  ... 
SIR10from4  SIR3from4  SIR4from5  SIR9from5  SIRlOfromS  SIR7from6  . 
SIRllfromG  SIR12from6  SIRllfrom7  SIR12from7  SIR13from7  ... 
SIR8from7  SIR12from8  SIR13from8  SIR14from8  SIR8from9  SIR13from9 


SIRl 

SIRl 

SIRl 

SIRl 

SIRl 

SIRl 

SIRl 


4from9  SIR15from9  S 
6fromll  SIR17fromll 
8froml2  SIR13froml2 
8froml4  SIR19froml4 
9froml5  SIR20froml5 
8froml7  SIR21froml8 
9from20  1  ]  )  .  '^  (-1)  ; 


IR9fromlO  SIR14fromlO  SIRl 


SIR12fromll 

SIR17froml3 

SIR20froml4 

SIR21froml6 

SIR18froml9 


SIR16froml2  S 
SIR18froml3  S 
SIR13froml4  S 
SIR17froml6  S 
SIR21froml9  S 


SfromlO  . . . 
IR17froml2  . . . 
IR19froml3  . . . 
IR14froml5  . . . 
IR21froml7  . . . 
IR21from20  . . . 


s  = 

[1 

112 

2  2 

2  3 

3  3  4 

4  4  4  5  5 

56667777888. 

9 

9  9  9 

10  1 

0  10 

11  11 

11  12  12 

12  12  13  13  13  14  14  14 

14  ... 

15 

15  15 

16 

16  17 

17  1 

8  19  19  20 

20  21]  ; 

D  = 

[2 

6  7  6 

7  8 

3  7 

8  9  8 

9  10  3  4 

9  10  7  11  12  11  12  13  8 

12  13  14 
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8  13  14  15  9  14  15  16  17  12  16  17  18  13  17  18  19  18  19  20  13 
14  19  20  21  17  21  18  21  18  21  21  19  20]; 

DG  =  sparse ( S , D, W) ; 

h  =  view (biograph (DG, [],' ShowWeights' on' ) ) 

[dist,  path,  pred]  =  graphshortestpath (DG,  1,  21)  ; 

dist 

path 
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